2009-08-13 13 views
0

longitude. J'ai un fichier .flt pour l'hémisphère nord qui contient toutes ces informations et je sais comment le lire, je ne sais pas comment faire correctement la projection de la carte. Je suis capable de faire la projection en utilisant LambertConformal (originLat, originLon, stdParallel1, stdParallel2) mais je ne sais pas comment utiliser les valeurs x, y projetées dans la grille. Les valeurs de grille sont Nx = 1073 et Ny = 689 et les valeurs projetées que j'obtiens pour lat = 33,92, lon = -84,33 sont x = 995 y = 1035. Je sais que les valeurs projetées sont correctes, mais je ne sais pas comment lire la grille avec eux.projection conique conforme Lambert pour obtenir les valeurs x, y d'une grille en java

Peut-être que quelqu'un peut m'aider à comprendre les valeurs suivantes, aussi.

:grid_mapping_name = "lambert_conformal_conic"; 
:standard_parallel = 25.0; // double 
:longitude_of_central_meridian = 265.0; // double 
:latitude_of_projection_origin = 25.0; // double 
:earth_shape = "Earth spherical with radius specified by producer"; 
:spherical_earth_radius_meters = 6371200.0; // double 
:GRIB_param_Dx = 5079.4062; // double 
:GRIB_param_Dy = 5079.4062; // double 
:GRIB_param_GDSkey = 1158554943; // int 
:GRIB_param_La1 = 20.192; // double 
:GRIB_param_LaD = 25.0; // double 
:GRIB_param_Latin1 = 25.0; // double 
:GRIB_param_Latin2 = 25.0; // double 
:GRIB_param_Lo1 = 238.446; // double 
:GRIB_param_LoV = 265.0; // double 
:GRIB_param_NpProj = "true"; 
:GRIB_param_Nx = 1073; // int 
:GRIB_param_Ny = 689; // int 
:GRIB_param_ProjFlag = 0; // int 
:GRIB_param_Quasi = "false"; 
:GRIB_param_ResCompFlag = 0; // int 
:GRIB_param_SpLat = -90.0; // double 
:GRIB_param_SpLon = 0.0; // double 
:GRIB_param_Winds = "True"; 
:GRIB_param_grid_name = "Lambert Conformal"; 
:GRIB_param_grid_radius_spherical_earth = 6371200.0; // double 
:GRIB_param_grid_shape = "Earth spherical with radius specified by producer"; 
:GRIB_param_grid_shape_code = 1; // int 
:GRIB_param_grid_type = 30; // int 
:GRIB_param_grid_units = "m"; 

Répondre

1

Avez-vous essayé de poster sur le site Unidata? Ils répondront probablement mieux à ce genre de questions. Dans le passé, j'ai reçu de très bonnes réponses de leur part. -Fred

+0

Cela semble être une suggestion raisonnable! Je vais devoir regarder dans leurs listes de diffusion, merci. – ravun