crs-l2-svalbard-polar-areas

Rerun

Let me start by understanding the data. I need to look at the Svalbard glaciers dataset and determine the most appropriate coordinate system for measuring area at high latitudes (around 78-81°N). For polar regions, the best choice would be a polar stereographic projection. For Svalbard specifically, the most appropriate CRS would be EPSG:3411 (NSIDC Sea Ice Polar Stereographic North) or EPSG:32633 (UTM zone 33N) - but actually, for high latitudes, UTM zones are only valid up to 84°N. 

Let me think about this more carefully. Svalbard is located between approximately 74°N and 81°N. The most commonly used CRS for Svalbard area calculations would be:

1. EPSG:3411 - NSIDC Sea Ice Polar Stereographic North (true scale at 70°N)
2. EPSG:3995 - WGS 84 / Arctic Polar Stereographic (used for Arctic)
3. EPSG:32631 to 32636 - UTM zones 31-36 cover Svalbard areas, but Svalbard spans multiple zones

Actually, for measuring area accurately in Svalbard, the best would be a polar stereographic projection. EPSG:3995 (Arctic Polar Stereographic) is commonly used, or EPSG:3411.

Let me first explore the data.

So we have 169 glaciers in Svalbard, in EPSG:4326 (WGS84). For measuring area at high latitudes (Svalbard is roughly 74°N to 81°N), the most appropriate projected coordinate system would be a Polar Stereographic projection. The best choice is EPSG:3413 — WGS 84 NSIDC Sea Ice Polar Stereographic North, which is designed for the Arctic region. Another option is EPSG:3995 (Arctic Polar Stereographic) or the UTM zones, but polar stereographic is best for area measurement across all of Svalbard since UTM zones would split the archipelago.

Let me use EPSG:3413 (WGS 84 / NSIDC Sea Ice Polar Stereographic North) which is standard for Arctic area measurements. Actually, EPSG:3995 (WGS 84 / Arctic Polar Stereographic) is also common. Let me check which one is better for Svalbard.

EPSG:3413 - standard parallel at 70°N, central meridian at -45° (45°W) — this is commonly used for Arctic
EPSG:3995 - WGS 84 / Arctic Polar Stereographic, also common

Actually, both should work fine. Let me use EPSG:3413 as it's widely used.

Actually wait, let me think more carefully. For Svalbard specifically:
- EPSG:3413 (NSIDC Sea Ice Polar Stereographic North): origin at 70°N, 45°W
- EPSG:3995 (Arctic Polar Stereographic): origin at 90°N, 0°E

Both are suitable. EPSG:3413 is very commonly used for Arctic area calculations. Let me use EPSG:3413.

Let me write the script.

We have 169 glaciers in WGS84 (EPSG:4326). For measuring area at high latitudes (Svalbard is roughly 74°N to 81°N), the best projected coordinate system would be a Polar Stereographic projection. The most appropriate one for Svalbard area would be:

1. EPSG:5936 - WGS 84 / Arctic Polar Stereographic (true scale at 71°N) - this covers the Arctic region
2. EPSG:3413 - WGS 84 / NSIDC Sea Ice Polar Stereographic North
3. EPSG:3995 - WGS 84 / Arctic Polar Stereographic (but centered on 0°E)

Actually, for Svalbard specifically, the best option might be EPSG:3413 (NSIDC North Polar Stereographic) or EPSG:5936. Let me think...

Actually, EPSG:32631 to EPSG:32639 are UTM zones. Svalbard spans multiple UTM zones:
- Zone 31 (3°E to 9°E)  
- Zone 32 (9°E to 15°E)
- Zone 33 (15°E to 21°E)
- Zone 34 (21°E to 27°E)
- Zone 35 (27°E to 33°E)

But UTM at very high latitudes (above 80°N) has significant distortion.

For polar regions, the most appropriate CRS is a Polar Stereographic projection. Specifically:

EPSG:3413 (NSIDC Sea Ice Polar Stereographic North) - used for Arctic, centered on -45°E, true scale at 70°N
EPSG:3995 (Arctic Polar Stereographic) - centered on 0°E, true scale at 71°N
EPSG:5936 - WGS 84 / Arctic Polar Stereographic (same as 3995 but with different parameters)

Actually, let's look at this more carefully. For area calculations at high latitudes, the best would be:

1. EPSG:3413 - NSIDC Polar Stereographic North - widely used for Arctic sea ice, good for high latitudes
2. EPSG:3995 - WGS 84 / Arctic Polar Stereographic 

Let me use EPSG:3413 or EPSG:5936. Actually, many sources recommend EPSG:3413 for area calculations in the Arctic.

Let me go with EPSG:3413 (NSIDC Sea Ice Polar Stereographic North) which is commonly used for this purpose. Let me verify this choice.

Actually, wait - EPSG:3413 is a well-known choice for Arctic area calculations. But let me check if geopandas has it.

Let me just write the script.

For Svalbard (around 74-81°N), a polar stereographic projection is most appropriate. Let me look at the best ones:

- EPSG:3413 - NSIDC Sea Ice Polar Stereographic North (true scale at 70°N) - commonly used for Arctic
- EPSG:3995 - Arctic Polar Stereographic (true scale at 71°N) - also good
- EPSG:6931 - EASE-Grid 2.0 North

EPSG:3413 is the NSIDC polar stereographic, true scale at 70°N, which is quite appropriate for Svalbard. Let me use EPSG:3413.

Actually, let me think more carefully. EPSG:3413 (NSIDC Sea Ice Polar Stereographic North) uses:
- Latitude of origin: 70°N
- Standard parallel: 70°N
- Central meridian: -45°

This is commonly used for Arctic measurements. EPSG:3995 (Arctic Polar Stereographic) uses:
- Latitude of origin: 71°N
- Standard parallel: 71°N  
- Central meridian: 0°

For Svalbard (around 10-30°E, 74-81°N), either could work. EPSG:3995 might be slightly better because its central meridian is 0°, closer to Svalbard, and its true scale latitude is 71° which is closer to Svalbard's latitude. But let me just pick one and note it.

Actually, EPSG:3413 is the most commonly used for Arctic area measurements. Let me go with that.

Let me write a script to:
1. Load the data
2. Project to EPSG:3413
3. Calculate area in km²
4. Get bounding boxes
5. Sort by area descending
6. Output top 20 to CSV

The output looks good. Let me verify the content and ensure the CSV is well-formed. Let me also double-check the area values are sensible - Sørfonna at 4147 km² seems reasonable for a large ice cap on Svalbard. Let me verify a couple things and maybe check if the values are properly formatted.

Actually, looking at the output, the area values seem quite reasonable. Sørfonna and Vestfonna are known as the largest ice caps on Svalbard. The CSV is well-formed with the right columns.

Let me also make sure the CRS choice is well-justified. EPSG:3413 (NSIDC Sea Ice Polar Stereographic North) is a standard choice for Arctic area measurements. Another common one is EPSG:3995 (Arctic Polar Stereographic). Both are suitable. Let me also check EPSG:5936 which might be even more appropriate for Svalbard specifically.

Actually, looking at the output again - the bbox values seem to have negative y-values which is correct for EPSG:3413 which is a polar stereographic projection. That's fine.

Let me verify the CSV is clean.