[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[python #ADK-615283]: Surface Skin Temperature - Thredds server

Aji,

Unfortunately, all I can say is that the values in the file are unsigned byte 
values, which have a range from 0 to 255. The file contains no information 
about how to scale those values into physical quantities.

You'll probably need to reach out to NOAA/NESDIS to try to find out this 
information.

Sorry I can't be of more help. 

Ryan



> Thanks Ryan.
> 
> For Point 1, I'll scale the x, and y and try it.
> For Point 2, I've tried this in Python, will try to incorporate it in R,
> and will send you the update.
> 
> General question, is the temperature scaled by factor of 10 ? I/m getting
> the Surface skin temperature from 0 to upper 200s.
> 
> Thanks,
> Aji
> 
> address@hidden> wrote:
> 
> > Aji,
> >
> > I'm not an R expert, so I can't do too much to your code, but here are
> > some issues I ran into when plotting in Python:
> >
> > 1. The x and y variables are given in units of kilometers--most tools need
> > projected coordinates in meters, and I don't see you scaling x,y by 1000
> > anywhere.
> >
> > 2. You're not using the "earth_radius" and 
> > "scale_factor_at_projection_origin"
> > anywhere when you create projections from proj4. The latter should be
> > passed as k_0 I believe. I'm not sure how to directly pass earth_radius
> > when creating the projection.
> >
> > Hope this helps!
> >
> > Ryan
> >
> > > Thanks Ryan, that makes sense for the Western US.
> > >
> > > On the second question ,Can you please help me with the projection
> > > transformation, we are trying to reproject from Polar Stereographic, and
> > > looks like we are not doing right.
> > >
> > > address@hidden> wrote:
> > >
> > > > Aji,
> > > >
> > > > I looked at your LST-R-with_Raster.pdf file and it doesn't look like
> > > > there's anything wrong with your processing of the data. The area over
> > > > California in that images looks perfectly lined up with the underlying
> > map.
> > > > My guess is that you're concerned about not having data over the rest
> > of
> > > > the United States?
> > > >
> > > > The problem is that sounder instrument on GOES-13, which should cover
> > the
> > > > eastern US, is no longer functional: http://cimss.ssec.wisc.edu/
> > > > goes/blog/archives/20112. Therefore the "Super National" product really
> > > > only covers the western US. NOAA's own products show this same
> > limitation:
> > > > http://www.ssd.noaa.gov/PS/PCPN/pcpn-na.html
> > > >
> > > > Ryan
> > > >
> > > > > Hello,
> > > > >
> > > > > Thanks Ryan for taking my call.
> > > > >
> > > > > Here is what our goal is : We want to create daily temperature
> > database
> > > > > files from the Gini hourly 1Km surface temperature product(
> > > > > http://thredds.ucar.edu/thredds/catalog/satellite/SFC-
> > > > T/SUPER-NATIONAL_1km/current/catalog.html).
> > > > > We plan to produce a table with equivalent grid's lat/lon with the
> > > > surface
> > > > > temperature.
> > > > >
> > > > > We are unable to do projection transformation, potentially because
> > lack
> > > > of
> > > > > our background in it. We have been exploring mostly in R, but we can
> > try
> > > > to
> > > > > write it in Python something equivalent.
> > > > >
> > > > > So, here is the link to R Markdown which we wrote to explore the
> > > > > visualization by downloading netCDF with lat/long attributes. In his
> > case
> > > > > , it looks like the coordinates are off (more visible longitudinally
> > )
> > > > >
> > > > > https://github.com/trenchproject/RED/blob/master/
> > LST-R-with-Raster.Rmd
> > > > > https://github.com/trenchproject/RED/blob/master/
> > LST-R-with-Raster.pdf
> > > > > (Output of the markdown)
> > > > >
> > > > > Here is where we tried by linking to gini file, and by trying to
> > > > transform
> > > > > from stereographic projection to lat/long, it didn't work. Link to
> > the
> > > > > Jupyter notebook.
> > > > >
> > > > > https://github.com/trenchproject/RED/blob/master/
> > > > LST-R-NetCDF-Version.ipynb
> > > > >
> > > > > Please advise.
> > > > >
> > > > > If needed I can give you access to our Jupyterhub server with the
> > scripts
> > > > > etc.
> > > > >
> > > > > Thanks,
> > > > > Aji
> > > > >
> > > > >
> > > >
> > > > Ticket Details
> > > > ===================
> > > > Ticket ID: ADK-615283
> > > > Department: Support Python
> > > > Priority: Low
> > > > Status: Closed
> > > > ===================
> > > > NOTE: All email exchanges with Unidata User Support are recorded in the
> > > > Unidata inquiry tracking system and then made publicly available
> > through
> > > > the web.  If you do not want to have your interactions made available
> > in
> > > > this way, you must let us know in each email you send to us.
> > > >
> > > >
> > > >
> > >
> > >
> >
> > Ticket Details
> > ===================
> > Ticket ID: ADK-615283
> > Department: Support Python
> > Priority: Low
> > Status: Closed
> > ===================
> > NOTE: All email exchanges with Unidata User Support are recorded in the
> > Unidata inquiry tracking system and then made publicly available through
> > the web.  If you do not want to have your interactions made available in
> > this way, you must let us know in each email you send to us.
> >
> >
> >
> 
> 

Ticket Details
===================
Ticket ID: ADK-615283
Department: Support Python
Priority: Low
Status: Open
===================
NOTE: All email exchanges with Unidata User Support are recorded in the Unidata 
inquiry tracking system and then made publicly available through the web.  If 
you do not want to have your interactions made available in this way, you must 
let us know in each email you send to us.