Suite

Conversion de la conique conforme de Lambert (LCC) en UTM

Conversion de la conique conforme de Lambert (LCC) en UTM


j'en ai deuxnumpymaillages (X, Y les deux 2D) dans la conique conforme de Lambert avec les informations suivantes :

J'aimerais les projeter directement (ou indirectement si ce n'est pas possible) vers UTM. j'ai les deuxpyprojetgdalinstallé, mais n'a pas été en mesure de comprendre comment les utiliser. J'ai trouvé des solutions pour la conversion de LCC en WGS84 ici et ici, mais je n'ai pu utiliser aucune d'entre elles car je n'ai pas toutes les informations nécessaires.

Même si d'une manière ou d'une autre je peux convertir deLCCàWGS84, alors je peux utiliser PythonUTMpaquet etutm.from_latlonpour obtenirUTMen dehors de ça.

J'ai trouvé la ligne de code suivante d'ici :

awips221 = Proj(proj='lcc', R=6371200, lat_1=50, lat_2=50,lon_0=-107, ellps="clrk66")

Et puis je peux utilisertransformerdepyproj. Dans ce cas, estlat_1etlat_2tous les deux25.0dans mon cas? Que seraitlon_0etellps?

Utilisantgdalinfonous pouvons voirlCCetWGS84projections pour le centre et les quatre coins de la grille dans l'image ci-dessous :


Je ne sais pas ce que la NOAA pense être les bonnes coordonnées, mais je n'ai aucun problème à charger le fichier dans QGIS ou à le reprojeter sur WGS84 avec gdalwarp :

QGIS utilise cette chaîne de projection personnalisée :

+proj=lcc +lat_1=25 +lat_2=25 +lat_0=25 +lon_0=265 +x_0=0 +y_0=0 +a=6371229 +b=6371229 +units=m +no_defs

Vous pouvez utiliser la même chaîne avec pyproj.Proj().

où vous pouvez remplacer en toute sécuritélon_0=265aveclon_0=-95.

Notez que l'étendue en degrés n'est pas la même que pour reprojeter les coordonnées du coin, car l'image est fortement courbée.

Je ne sais pas pourquoi vous voulez transformer les données en UTM. Il traverse plusieurs zones UTM (7N à 21N). Vous pouvez prendre celui du milieu (15N) et subir des distorsions sur les côtés, ou couper l'image le long des limites de la zone UTM et effectuer des calculs dans chaque zone UTM séparément.


J'ai essayé cette chaîne dans pyproj :

Proj(proj='lcc', R=6371200, lat_1=25, lat_2=25,lon_0=-95, ellps="clrk66")

et mes résultats étaient étranges. quand j'ai changé les noms de champs lat:

Proj(proj='lcc', R=6371200, lat_0=25, lat_1=25,lon_0=-95, ellps="clrk66")

J'ai obtenu les résultats que j'attendais.


Voir la vidéo: Map Projections Part 2: Conic Projections