gis - Reading GlobSnow HDF4 files in R -


i trying read .hdr files in r. files come globsnow , can found @ http://www.globsnow.info/swe/archive_v2.0/2013/l3a_daily_swe_hdf/

i have used gdalutils package, seems work modis data, not working globsnow data. have tried this:

sd <- get_subdatasets("globsnow_swe_l3a_20130101_v2.0.hdf") gdal_translate(sd[2], sds=t, dst_dataset = "hdftrial.tif")

however returns null.

i have tried installing gdal through homebrew, unsure how make work in r. opened data in hdfview , data seems right. there doesn't seem easy way access data in r. hoping raster or matrix .hdf file.

the null comes out me too, follow raster("hdftrial.tif") , loads raster in no problem. null comes because have no subdatasets (sds), intentionally calling single dataset, appears hdf not read in subdatasets if not declare.

library(gdalutils) library(raster) sd <- get_subdatasets("c:/_localdata/globsnow_swe_l3a_20130101_v2.0.hdf") gdal_translate(sd[2], dst_dataset = "hdftrial.tif")  plot(raster("hdftrial.tif"))

when use gdal_translate pull each individual dataset both there, throws null if there not "subdatasets".

from gdal_translate:

 if (output_raster) {         if (missing(sds)) {             return(brick(dst_dataset))         }         else {             if (!sds) {                 return(brick(dst_dataset))             }             else {                 return(null)             }         }     }     else (return(null))

update:

gdal_translate("c:/_localdata/globsnow_swe_l3a_20130101_v2.0.hdf", sds=t, of="gtiff", dst_dataset = "hdftrial.tif")

you not require get_subdatasets, internally. declaring sds=t create individual file each subdataset. still throws null, trying decipher why. file _#.tif each subdataset exists within dataset.


Comments

Post a Comment

Popular posts from this blog

javascript - Create a stacked percentage column -

i have collection of data : var data = [ {"p301a":"10","p301b":"7","p301c":"7","p301d":"3","p301e":"8","p301f":"8","p301g":"8","p301h":"8","p301i":"8","p301j":"8","p301k":"8","p301l":"8","p301m":"8","p301n":"8","p301o":"8","age":"31-40 years","profesion":"2","position":"2"}, {"p301a":"5","p301b":"4","p301c":"4","p301d":"4","p301e":"4","p301f":"4","p301g":"4","p301h":"4","p301i":"4","p301j":"4","p301k":"4...
Read more

Optimising Firebase database by automatically overwriting data -

in firebase app have following structure: posts latest all users can write posts > latest . want limit entries 20 posts save space in database. , posts beyond 20 unnecessary never shown on homepage. how can limit posts 20, when user writes posts > latest posts @ end drop off (be deleted) automatically? this sample should you. you need : const max_log_count = 20; exports.removeold = functions.database.ref('/posts/latest/{postid}').oncreate(event => { const parentref = event.data.ref.parent; return parentref.once('value').then(snapshot => { if (snapshot.numchildren() >= max_log_count) { let childcount = 0; const updates = {}; snapshot.foreach(function(child) { if (++childcount <= snapshot.numchildren() - max_log_count) { updates[child.key] = null; } }); // update parent. removes children. ...
Read more

javascript - Angular UI-Grid customTemplate directive causing rows to load slowly/? -

Image
each grid cell can have more 1 customtemplate loaded. have directive holds watcher on model_col_field listen cell updates on scroll when dom rows updated , recompile celltemplate . without watcher , $compile rows repeat on scroll, showing same set of rows. is there more angular way of approaching or way improve performance? relevant controller code: create gridoptions.data & gridoptions.columndefs : var customtemp = '<div ng-repeat="opt in grid.appscope.colvalueoptionsmodel" ng-init="tempstring=model_col_field.clazz+opt.key;"><comparison-directive dbo="model_col_field" attobj="opt" fff="grid.appscope.getattributetemplate(tempstring)" /></div>'; $scope.gridoptions.columndefs = [ { name:'yindex',pinnedleft:true,celltemplate: customtemp, displayname:'index'}, { name:'col1',celltemplate: customtemp}, { name:'col2',celltemplate: c...
Read more