R2 Landbrugsareal (Agricultural Land) Data Catalog

Field boundaries, crop data, land use, and agricultural production data.

Frontend Metrics Supported

Available Datasets

Silver Layer

FVM Marker - Field Boundaries (617K rows/year)

Path: r2://landbruget-data/silver/fvmmarker{year}/*/data.parquet Years Available: 2008-2025

Schema (introspected):

field_id: string
area_ha: double
cvr_number: string
crop_code: int64
crop_name: string
block_id: string
geometry: binary (WKB)
year: int64
field_uuid: string
municipality: string
is_organic: bool
organic_conversion_date: date32
[20 columns total]

Agricultural Blocks

Path: r2://landbruget-data/silver/agriculturalblocks{year}/*/data.parquet

Cadastral (2.16M rows)

Path: r2://landbruget-data/silver/cadastral/*/data.parquet

BBR Buildings (579K rows)

Path: r2://landbruget-data/silver/bbrbuildings/*/joinedbuildings.parquet

DAGI Kommuner (99 rows)

Path: r2://landbruget-data/silver/dagi_kommuner/*/data.parquet

Gold Layer

Field Production (617K rows/year)

Path: r2://landbruget-data/gold/fieldproduction{year}/*/data.parquet

Schema (introspected):

field_id: string
block_id: string
cvr_number: string
year: int64
area_ha: double
crop_type: string
organic_farming: bool
landsdel_code: string
landsdel_name: string
dst_regions: string
kommune_name: string
yield_estimate_hkg_ha: double
production_estimate_hkg: double
field_uuid: string
primary_field_id: string
[18 columns total]

Common Queries

Get All Fields for a CVR

import duckdb
from common.storage.filesystem import setup_duckdb_cloud_auth

conn = duckdb.connect()
setup_duckdb_cloud_auth(conn)

# Read FVM marker data for 2024
df = conn.execute("""
    SELECT * FROM read_parquet('r2://landbruget-data/silver/fvm_marker_2024/2025-01-10/data.parquet')
""").df()

# Filter by CVR
company_fields = df[df['cvr_number'] == '31373077']
print(f"Total fields: {len(company_fields)}")
print(f"Total area: {company_fields['area_ha'].sum():.2f} ha")

Calculate Organic Farming Percentage by Municipality

# Group by municipality and organic status
organic_stats = df.groupby(['municipality', 'is_organic']).agg({
    'area_ha': 'sum'
}).unstack(fill_value=0)

organic_stats.columns = ['conventional', 'organic']
organic_stats['organic_pct'] = (
    organic_stats['organic'] /
    (organic_stats['conventional'] + organic_stats['organic']) * 100
)
organic_stats = organic_stats.sort_values('organic_pct', ascending=False)

Get Production Estimates for a Crop

# Read gold production data
production = conn.execute("""
    SELECT * FROM read_parquet('r2://landbruget-data/gold/field_production_2024/2025-01-10/data.parquet')
""").df()

# Filter for wheat
wheat = production[production['crop_type'].str.contains('hvede', case=False)]
print(f"Total wheat area: {wheat['area_ha'].sum():.2f} ha")
print(f"Average yield: {wheat['yield_estimate_hkg_ha'].mean():.2f} hkg/ha")
print(f"Total production: {wheat['production_estimate_hkg'].sum():.2f} hkg")

Crop Diversity Analysis

# Count unique crops per CVR
crop_diversity = df.groupby('cvr_number').agg({
    'crop_name': 'nunique',
    'area_ha': 'sum'
}).reset_index()
crop_diversity.columns = ['cvr_number', 'crop_count', 'total_area']
crop_diversity = crop_diversity.sort_values('crop_count', ascending=False)

Read Geometry with DuckDB Spatial

# Read with geometry via DuckDB
gdf = conn.execute("""
    SELECT * FROM read_parquet('r2://landbruget-data/silver/fvm_marker_2024/2025-01-10/data.parquet')
    WHERE cvr_number = '31373077'
""").df()

Join Keys

Data Quality Notes

FVM Marker

• Update frequency: Annual (new year data released ~Q2)
• Coverage: All registered agricultural fields in Denmark
• Years available: 2008-2025 (varies by data source)
• Geometry: WKB format, EPSG:4326

Field Production

• Update frequency: Annual after harvest statistics
• Coverage: Production estimates for all registered fields
• Caveat: Yield estimates are modeled, not measured

Cadastral

• Update frequency: Weekly
• Coverage: All cadastral parcels in Denmark
• Note: Use for property boundaries, not field boundaries

Related Skills

• okonomi/ - Subsidies linked to fields via CVR and markbloknummer
• miljo/ - Environmental data at field level (pesticides, nitrogen)
• husdyr/ - Livestock density per land area

R2 Paths Reference

# List available FVM marker years
rclone lsd r2:landbruget-data/silver/ | grep fvm_marker

# List field production years
rclone lsd r2:landbruget-data/gold/ | grep field_production

# List cadastral snapshots
rclone lsd r2:landbruget-data/silver/cadastral/

# List DAGI municipality data
rclone lsd r2:landbruget-data/silver/dagi_kommuner/

Yearly Data Pattern

For datasets with yearly data, use this pattern:

# Query specific year
year = 2024
path = f'silver/fvm_marker_{year}/*/data.parquet'

# Query all years (use with caution - large data)
path = 'silver/fvm_marker_*/*/data.parquet'