# YellowSea - 3413 Doesn't work correctly
# fileName = 'ASA_IMM_1PNPDE20080202_021339_000000812065_00361_30976_3972.N1'
# fileName = 'ASA_WSM_1PNPDK20120320_135328_000000923112_00427_52596_1286.N1'

# Baikal - same thing
# fileName = 'ASA_WSM_1PNPDE20100912_025847_000000922092_00476_44618_1878.N1'

# Kuroshio - also cropping is not correct enough
# fileName = 'ASA_WSM_1PNPDK20120108_125143_000000923110_00254_51561_3395.N1'

# Laptev - awful cropping in 3413 around 150E
# (134.73301433365424, 158.69928524693626) (67.91966824291822, 78.468371541120206)

# Looks like the problem is with lons around 110-140 East
# http://epsg.io/3413 notes that the EPSG3413 is a Cartesian 2D CS for north polar azimuthal lonO 45°W.
# Axes: X,Y. Orientations: X along 45°E, Y along 135°E meridians. UoM: m.
# That could be the case!!!! Y along 135°E meridians!!!!!!
# In any case this coordinate system is used for above 60°N only!!!!
# YES looks like that this was the case, for 3413 lat_0=90 and lon_0=-45
# SO THE PROJ DOESN'T SEE some parts of the image
# If you specify lat_0/lon_0 for the center of the image everything works great
#         lat_0 = (up + down) / 2
#         lon_0 = (right + left) / 2
# But that wouldn't be 3413 anymore :)