From:	Sergio Andrés Marchant Rojas <>
To:	Jason Roberts <>
Cc:	Sergio Marchant <>,
Subject:	Re: [mget-help] How to create Hovmoller plots and other figures
Date:	Sun, 4 Jul 2010 00:21:47 -0400

Dear Jason et al. and MGET users,

Thanks for the information. Now I know where to look for tools to make
those graphs. I have another question that I could use some help too.
I am trying to import AVISO MSLA merged files but am having troubles
to find the information that is required to use the tool 'Convert 2D
Variable in NetCDF to ArcGIS Raster'. I am specifically confused
regarding the x coordinate, y coordinate, cell size and projection
conversions fields.

I tried to find those variables in the NetCDF header however that was
confusing as also when trying to understand the section 2.1.7.2. ‘Maps
projection’ of the SSALTO/DUACS User Handbook.

Here are the variables that I used (my best guess), however I am not
sure why it does not correctly align with my SST dataset imported from
Pathfinder v5.0 (WGS_84 projection)

Variable name: Grid_0001
x coordinate: -180
y coordinate: -90
cell size: 0.333333333333333

Preconversion:

Transpose and flip

Postconversion:

Define coordinate system: Mercator(Sphere)
Project to new coordinate system: GCS_WGS_1984
Projection resampling: NEAREST
cell size for projected coord system: 0.333333333333333

The final map has the lowest value of -123.729 and the highest of
1.84467e+019. This values should be in cm however and I think they are
not correct. Thanks again for your help.

Here is a copy of the NetCDF header:

netcdf F:\OCEANO
DATA\GeoEcoTemp\tmpfynjut\inputFile\dt_upd_global_merged_msla_h_19921014_19921014_20100503
{
dimensions:
        NbLatitudes = 915 ;
        NbLongitudes = 1080 ;
        GridDepth = 1 ;
        LatLon = 2 ;
variables:
        int LatLon(LatLon) ;
                LatLon:_FillValue = 2147483647 ;
                LatLon:long_name = "No sense but necessary for some automatic 
tools" ;
                LatLon:units = "count" ;
        double NbLatitudes(NbLatitudes) ;
                NbLatitudes:_FillValue = 1.84467440737096e+019 ;
                NbLatitudes:long_name = "Latitudes" ;
                NbLatitudes:units = "degrees_north" ;
        double NbLongitudes(NbLongitudes) ;
                NbLongitudes:_FillValue = 1.84467440737096e+019 ;
                NbLongitudes:long_name = "Longitudes" ;
                NbLongitudes:units = "degrees_east" ;
        double LatLonMin(LatLon) ;
                LatLonMin:_FillValue = 1.84467440737096e+019 ;
                LatLonMin:long_name = "Latitude/Longitude of south/ouest 
corner" ;
                LatLonMin:units = "degree" ;
        double LatLonStep(LatLon) ;
                LatLonStep:_FillValue = 1.84467440737096e+019 ;
                LatLonStep:long_name = "latitude/longitude steps" ;
                LatLonStep:units = "degree" ;
        float Grid_0001(NbLongitudes, NbLatitudes) ;
                Grid_0001:_FillValue = 1.844674e+019f ;
                Grid_0001:long_name = "H" ;
                Grid_0001:units = "cm" ;
                Grid_0001:Date_CNES_JD = 15627. ;
                Grid_0001:date = "1992-10-14 00:00:00.000000 UTC" ;

// global attributes:
                :FileType = "GRID_DOTS_MERCATOR" ;
                :OriginalName = 
"dt_ref_global_merged_msla_h_19921014_19921014_20100503.nc" ;
                :CreatedBy = "SSALTO/DUACS" ;
                :CreatedOn = "05-MAY-2010 22:30:13:000000" ;
                :title = "SSALTO/DUACS - DT MSLA - Merged Product - 
Homogeneous
Global Processing" ;
}

Regards,

Sergio.


2010/7/1, Jason Roberts 
<>:
> Hi Sergio,
>
>
>
> Thanks for contacting us about this. Unfortunately I don't know of an easy
> way to produce these kinds of plots using ArcGIS. When I've done similar
> plots in the past I've used R or MATLAB. I don't have any example code handy
> but you could probably find some just by looking around. If you prefer to
> use Python, I suggest reading the data using GDAL's Python bindings (or an
> internal function in MGET that wraps them), performing any needed
> statistical summarization using standard numpy array manipulation, and then
> using the matplotlib package to create the plots. I'm sorry I can't be more
> help than that.
>
>
>
> We do have a long term interest in being able to produce these kinds of
> plots with MGET. We are currently rewriting MGET's raster processing modules
> to be aware of 3D and 4D data. Once that is done, we will build basic tools
> that allow conversion, sampling, and building of climatologies for 3D and 4D
> data. Having tools to produce Hovmoller plots would also be useful. I will
> add that to our "to do" list.
>
>
>
> Best regards,
>
>
>
> Jason
>
>
>
> From: Sergio Marchant 
> [mailto:]
> Sent: Thursday, July 01, 2010 11:28 AM
> To: 
> 
> Subject: [mget-help] How to create Hovmoller plots and other figures
>
>
>
> Dear Dr. Roberts,
>
>
>
> I am a Biological Science graduate student, and I am new using ArcGIS and
> the Marine Geospatial Ecology tools, however I am progressing quickly thanks
> to the excellent supporting documentation and examples. I would like to know
> if it is possible to create a shaded Hovmoller plot that describes monthly
> average values in the x axis (time) and other parameters (i.e Sea Heigh,
> Winds, currents, Chlorophyl, Sea surface temperature) in the Y axis . I
> already download monthly average data that I need in hdf and ncdf from
> several satellites and I thanks to MGET I can import the data in ArcGIS,
> however I am having trouble in figure out how to display all this
> information (temporal trends) in a scientific summarized way. Here are some
> examples of figures that I would like to include in my research.
>
>
>
>
>
> Thanks in advance.
>
>
>
> Sergio Marchant
>
> 
>
>

< date >
< thread >

Archives powered by MHonArc.