It is currently Mon, 18-11-19, 9:20 GMT

All times are UTC




Post new topic Reply to topic  [ 148 posts ]  Go to page Previous  1, 2, 3, 4, 5, 6 ... 10  Next
Author Message
 Post subject:
PostPosted: Sun, 28-10-07, 22:28 GMT 
Offline
User avatar

Joined: Thu, 25-10-07, 15:20 GMT
Posts: 992
Location: NE PA, USA
I am also seeing a LOT of areas that should not have water, at least there is nothing on the base texture that indicates there should be water. Just east of the Andes, there are a lot of areas. And that's not dessert. This arcexplorer might come in handy. I think that before any attempt to process the swbd files into a spec map, areas that should not have water should be identified and compared to the swbd file from that area to see if it also has the water.
cartrite


Top
 Profile  
 
 Post subject:
PostPosted: Mon, 29-10-07, 17:47 GMT 
Offline
User avatar

Joined: Thu, 25-10-07, 15:20 GMT
Posts: 992
Location: NE PA, USA
This first sample I isolated seems to suggest that these areas that are on the watermask file that should not be there are not on the swbd files. This sample is on tile tx_19_17 level5. I grayscaled the tile and enhanced a little to bring out the waterways better. The arrow points to the bad data.

Image

The spec map from the alpha channel.

Image

A screen capture from Arcexplorer with tiles covering 6s to 12s and 68w to 73w.

Image
The spec map from the alpha channel seems to have parts of a tributary not in the swbd? Not sure about all the white spots to the lower right of this last image.
cartrite


Top
 Profile  
 
 Post subject:
PostPosted: Mon, 29-10-07, 20:37 GMT 
Offline
Site Admin
User avatar

Joined: Fri, 31-08-07, 7:01 GMT
Posts: 4660
Location: Hamburg, Germany
Hi all,

after several days I now got a (quite enthusiastic) answer from Dr. Reto Stoeckli, lead scientist of the NASA BMNG project. Let me challenge first of all your German (he is Swiss). If you need help, no problem:

---------------------------------------------------------------
Fridger

Ihr macht sehr tolle Arbeit aus meine Bildern. Du solltest Zugriff zum 30m Landsat Inventory kriegen ... aber da habe ich auch nicht Zugriff, und Google Earth ist keine non-profit Organisation ;-)

Zu deinen texture - Tools sollte ich wahrscheinlich noch einen subsetting-server bauen damit man nicht die gesamten Daten runterladen muss ;-)

Hier noch einige Antworten:
>
> ----------------------------------------------
> Habt ihr einen aktzeptablen >= 86400x43200 /binaeren Watermask/ Datensatz irgenwie verfuegbar? Waere super...
> ----------------------------------------------

Leider habe ich keine anderen Watermasks (also es gibt eine oceanmask und eine watermask, das hast du gesehen?), und diejenigen die ich bieten kann sind aus den MODIS land team Datenbanken rausgelesen. Die Maske ist teilweise sehr fragwuerdig. Die Daten in der Arktis (>80N) sind noch schlechter da ich dort auf die Landklassifizierung angewiesen war (landcover = 0).

Kannst du solche Masks aus SRTM machen? Ich waere einen Abnehmer!!!

> Eine kurze weitere Frage:
>
> Gibt es bei Euch (BMNG) noch richtig grosse Packeis-layers fuer die Polarregionen? Sie existierten in geringerer Aufloesung in den alten BlueMarble TIF Originalen als separate Layers. Die BMNG textures sind alle ohne Packeis publiziert! Gibt's da einen hires Ausweg?
>

Es gibt einige fast eingefrorene (wegen Zeitmangels meinerseits) BMNG Projekte:
- monatliche Eismaske + Eistemperatur
- monatliche Ocean Chlorophyl Gehalt
- Wolkenbilder (Feature Oriented, z.B. mit Hurricanes, oder Polarstuerme)

Jaja, du siehst! Pick dir einees aus und ich motiviere dich daran zu arbeiten! Die MODIS Daten sind alle gratis.

Keep me posted and keep up the fun stuff!
Reto



--
Dr. Reto Stockli

Atmospheric Sciences - Colorado State University
NASA Earth Observatory - Goddard Space Flight Center
Atmospheric and Climate Science - ETH Zurich

http://www.iac.ethz.ch/people/stockli
http://bluemarble.nasa.gov
http://www.yawah.com/bmng

---------------------------------------------------

Bye Fridger

_________________
Image


Top
 Profile  
 
 Post subject:
PostPosted: Mon, 29-10-07, 20:52 GMT 
Offline
Site Admin
User avatar

Joined: Thu, 30-08-07, 22:52 GMT
Posts: 2727
Location: France, South, not far from Montpellier
t00fri wrote:
Hi all,

after several days I now got a (quite enthusiastic) answer from Dr. Reto Stoeckli, lead scientist of the NASA BMNG project. Let me challenge first of all your German (he is Swiss). If you need help, no problem:...


Fridger, BabelFish translation is a bit chaotic... :? Could you summarize please? :oops:


Top
 Profile  
 
 Post subject:
PostPosted: Mon, 29-10-07, 21:32 GMT 
Offline
User avatar

Joined: Thu, 25-10-07, 15:20 GMT
Posts: 992
Location: NE PA, USA
Quote:
Kannst du solche Masks aus SRTM machen? Ich waere einen Abnehmer!!!
Quote:
Can you make such MASKs from SRTM? I would be a customer!!!
:wink:
If the watermask used the MODIS landteam database then the swbd should be okay?
cartrite


Top
 Profile  
 
 Post subject:
PostPosted: Mon, 29-10-07, 22:32 GMT 
Offline
Site Admin
User avatar

Joined: Fri, 31-08-07, 7:01 GMT
Posts: 4660
Location: Hamburg, Germany
ElChristou wrote:
t00fri wrote:
Hi all,

after several days I now got a (quite enthusiastic) answer from Dr. Reto Stoeckli, lead scientist of the NASA BMNG project. Let me challenge first of all your German (he is Swiss). If you need help, no problem:...


Fridger, BabelFish translation is a bit chaotic... :? Could you summarize please? :oops:


C'est trop, mon ami ...
Alors, il faut apprendre l'Allemand (juste un peu) ;-)

F.

_________________
Image


Top
 Profile  
 
 Post subject:
PostPosted: Mon, 29-10-07, 22:34 GMT 
Offline
Site Admin
User avatar

Joined: Fri, 31-08-07, 7:01 GMT
Posts: 4660
Location: Hamburg, Germany
There was a second mail coming in that is way more interesting:

----------------------------------------------------------
Fridger,

Uebrigens, benuetzt du BMNG version 2 mit der enhanced Antarktis?

ftp://bluemarble.nascom.nasa.gov/bmng/

Falls du nicht einloggen kannst, musst du mir die IP Adresse deines Computers geben damit wir dich freischalten koennen.

Reto

-------------------english-------------------------
Fridger,

actually, do you use BMNG version 2 with enhanced Antarctica?

ftp://bluemarble.nascom.nasa.gov/bmng/

If you can't login, you got to send me your IP address so we can enable your access

Reto
-----------------------------------------------------

I cannot yet get in...

Bye Fridger

_________________
Image


Last edited by t00fri on Mon, 29-10-07, 22:49 GMT, edited 2 times in total.

Top
 Profile  
 
 Post subject:
PostPosted: Mon, 29-10-07, 22:36 GMT 
Offline
Site Admin
User avatar

Joined: Fri, 31-08-07, 7:01 GMT
Posts: 4660
Location: Hamburg, Germany
cartrite wrote:
Quote:
Kannst du solche Masks aus SRTM machen? Ich waere einen Abnehmer!!!
Quote:
Can you make such MASKs from SRTM? I would be a customer!!!
:wink:
If the watermask used the MODIS landteam database then the swbd should be okay?
cartrite


Perfect, cartrite!

I am still at a loss about the conclusions...

Bye Fridger

_________________
Image


Top
 Profile  
 
 Post subject:
PostPosted: Mon, 29-10-07, 22:42 GMT 
Offline
Site Admin
User avatar

Joined: Fri, 31-08-07, 7:01 GMT
Posts: 4660
Location: Hamburg, Germany
Did everyone get this?

I was asking Reto about the existence of polar ice masks, as they were shipped in lower resolution with the old BlueMarble data. His answer

------------------original, german --------------
Es gibt einige fast eingefrorene (wegen Zeitmangels meinerseits) BMNG Projekte:
- monatliche Eismaske + Eistemperatur
- monatliche Ocean Chlorophyl Gehalt
- Wolkenbilder (Feature Oriented, z.B. mit Hurricanes, oder Polarstuerme)
----------------------------------------------------

------------------english------------------------
There are some almost frozen BMNG projects (because of my lack of time):
  • monthly ice mask + ice temperature
  • monthly chlorophyl content in oceans
  • cloud images (feature oriented, e.g. with hurricans or polar storms!)

--------------------------------------------------

Bye Fridger

_________________
Image


Top
 Profile  
 
 Post subject:
PostPosted: Mon, 29-10-07, 22:55 GMT 
Offline
User avatar

Joined: Thu, 25-10-07, 15:20 GMT
Posts: 992
Location: NE PA, USA
Quote:
Perfect, cartrite!

Thanks to Google!!! I can make out some words though. I kind of understood those sentences.
cartrite


Top
 Profile  
 
 Post subject:
PostPosted: Tue, 30-10-07, 2:41 GMT 
Offline
User avatar

Joined: Sat, 13-10-07, 18:18 GMT
Posts: 373
English, German or French... It's all GREEK to me!!!

Sorry, couldn't resist. I'm going to TRASH my Vista installation this week I
believe. So, it may be a while until I can get organized again. I do have your
tools though, Dr. Schrempp, and I WILL learn how to use them as soon as I
get finished trashing this junk O.S.

Thanks, Brain-Dead

_________________
Brain-Dead Geezer Bob
Windows 8.1 Home Premium
AMD A8-5500 APU with Radeon HD Graphics
8.00 GB RAM, 64-bit Operating System
HP Hardware


Top
 Profile  
 
 Post subject:
PostPosted: Wed, 31-10-07, 21:59 GMT 
Offline
User avatar

Joined: Thu, 25-10-07, 15:20 GMT
Posts: 992
Location: NE PA, USA
Fridger,
I'm not sure yet but there seems to be some utilities called gdal that I'm going to try that may process the swbd files. There is a program called gdal_rasterize that seems to do just what we want. It sounds like if I create an ocean only specmap from the srtm84k file and use gdal_rasterize to "burn" the polygons from the swbd shape files onto the specmap. That is if I can get gdal to read the specmap. The program gdal_translate may do that. I downloaded the gdal utiliies from http://trac.osgeo.org/gdal/wiki/DownloadSource. I actually downloaded this to use on some files I create with ISIS but it seems to have the capability to process the swbd files as well. Fingers crossed.
cartrite


Top
 Profile  
 
 Post subject:
PostPosted: Thu, 01-11-07, 10:33 GMT 
Offline
Site Admin
User avatar

Joined: Fri, 31-08-07, 7:01 GMT
Posts: 4660
Location: Hamburg, Germany
Wow, cartrite,
that sounds interesting. I'll have a look tonight, too.

Bye Fridger

_________________
Image


Top
 Profile  
 
 Post subject:
PostPosted: Thu, 01-11-07, 11:42 GMT 
Offline
User avatar

Joined: Thu, 25-10-07, 15:20 GMT
Posts: 992
Location: NE PA, USA
This may take a while to ingest. In case you didn't see this already, here is a page with links that describe how gdal works. http://www.gdal.org/pages.html
EDIT:These same pages are also in the man folder from the source code.
EDIT: You will probably be able to understand those pages a lot better than I can since they are all about c code. I think the ocean only specmap needs to be read by gdal so all the shapefiles from the swbd can be scaled down and placed in the proper locations. This may be as simple as adding a header in the ocean-only specmap that gdal can read. Simpler yet, it may only require the keyword for datum to = WGS84. This is the info I think it needs? Maybe the .h files in the src code has the values for the field names that gdal looks for in the .prf file or image header.
Code:
xoffset=-43199
yoffset=-21599
pixelres=500 <meters>
datum=WGS84/EGM96? datum=WGS84
minlogitude=-180
maxlongitude=180
minlatitude=-90
maxlatitude=90
centerlongitude=0
cenerlatitude=0
blksize=?
samples=86400
lines=43200
projection=PlateCarree?
I think the program gdalwarp ( via .prf file ) can be used to enter this info into the ocean specmap? Then swbd files should have all the info needed in the .shp file to let gdal_rastereize scale down and place the polygons in the proper place "burned" to 255. I haven't tried anything yet because I'm not completely sure how to get gdal read both the specmap and swbd files as having the same coordinate system.

cartrite


Top
 Profile  
 
 Post subject:
PostPosted: Mon, 05-11-07, 19:24 GMT 
Offline
User avatar

Joined: Thu, 25-10-07, 15:20 GMT
Posts: 992
Location: NE PA, USA
I'm getting closer. I created a file called test.png with the gimp 3601x3601.
Then I ran gdal_translate with the following.
Code:
gdal_translate  -a_ullr -75 39 -74 40 test.png test.tif

Then ran gdal_rasterize with the following
Code:
gdal_rasterize -b 1 -burn 255 -l w075n39n -a WB w075n39n.shp test.tif

Which produced this.
This image was rgb but is full resolution. The reason the water is red is because I only burned 1 band. Should have burned 3.
Image
Scaled down to fit forum.
Image

For my second test, I created a black image at 21600x21600 for tile B1. Then I ran
Code:
cartrite@linux-o863:/srv/disk4/spectest> gdal_translate  -a_ullr -90 90 0 0 test1.png test1.tif
Input file size is 21600, 21600
0...10...20...30...40...50...60...70...80...90...100 - done.
cartrite@linux-o863:/srv/disk4/spectest> gdal_rasterize -b 1 -burn 255 -l w075n39n -a WB w075n39n.shp test1.tif
0...10...21...32...40...51...62...70...81...91...100 - done.

It "looks like it was placed in the right space. Can't show the full tile but this would be the full resolution specmap for that tile at 84k.
Image

cartrite


Top
 Profile  
 
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 148 posts ]  Go to page Previous  1, 2, 3, 4, 5, 6 ... 10  Next

All times are UTC


Who is online

Users browsing this forum: No registered users and 1 guest


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
cron
Powered by phpBB® Forum Software © phpBB Group