From:	"Jason Roberts" <>
To:	'Sergio Andrés Marchant Rojas' <>
Cc:	"'Sergio Marchant'" <>, <>
Subject:	RE: [mget-help] How to create Hovmoller plots and other figures
Date:	Mon, 12 Jul 2010 09:14:48 -0400

Hi Sergio,

Sorry for the slow response. I was on vacation next week.

I know the problem you are facing with trying to determine the map
projection parameters from the SSALTO/DUACs User Handbook. I carried out a
conversation with Aviso User services for almost a year trying to get to the
bottom of that problem without a satisfactory result. I have an answer that
should work for most people, but because the Aviso data are used to measure
global sea level rise, I wanted to try to be as exact as possible. I'll skip
the details for now, but let me know if you're interested.

The Aviso data are in a Mercator coordinate system. Therefore the corner
coordinates and cell size must be in linear units (meters), rather than
degrees. But before I go into these and the other parameters, may I ask
whether it is feasible for you to use the MGET tools that download Aviso
data directly over OPeNDAP, such as Download Aviso SSH Dataset to ArcGIS
Rasters? These tools are much easier to use than downloading and converting
the netCDFs yourself.

If you need to go the netCDF route, I can go ahead and tell you the
parameters. I would just do this anyway but I do not have them written
down--I have some code that calculates them--so it is a little bit of work
to do it.

Note that Aviso currently provides access to the DT version 2.1 data through
OPeNDAP but you can get DT 3.0 via the FTP files. That may be one reason why
you want the netCDF parameters. Aviso is planning to upgrade the OPeNDAP to
serve DT 3.0 but when this happens, the OPeNDAP will no longer allow
anonymous access. At that point, we will have to release a new version of
the MGET tool that supports authenticated OPeNDAP. You will have to enter
your Aviso user name and password into the tool. We are ready to do this but
are waiting for Aviso to provide us with details of how it should work.
Hopefully they will do that soon...

Jason

-----Original Message-----
From: Sergio Andrés Marchant Rojas 
[mailto:]
 
Sent: Sunday, July 04, 2010 12:22 AM
To: Jason Roberts
Cc: Sergio Marchant; 

Subject: Re: [mget-help] How to create Hovmoller plots and other figures

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_199
21014_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.