Extracting specific netcdf info and converting to GeoTIFF in python -


i attempting extract specific set of data netcdf file , convert said data geotiff.

so far have managed extract data want using netcdf4, data in file stored 1d arrays (lat, lon, data want) , assigning them 2d array. netcdf file working subsetted specific region. here @ loss.

i have slight understanding of how geotiff conversion works via have read @ these links:

https://borealperspectives.wordpress.com/2014/01/16/data-type-mapping-when-using-pythongdal-to-write-numpy-arrays-to-geotiff/

http://adventuresindevelopment.blogspot.co.uk/2008/12/create-geotiff-with-python-and-gdal.html

and here have currently:

import netcdf4 import numpy np osgeo import gdal osgeo import osr  #reading in data files , extracting said data ncfile = netcdf4.dataset("data.nc", 'r')  dataw = ncfile.variables["dataw"][:] lat = ncfile.variables["latitude"][:] long = ncfile.variables["longitude"][:]   n = len(dataw) x = np.zeros((n,3), float)  x[:,0] = long[:] x[:,1] = lat[:] x[:,2] = dataw[:]  nx = len(long) ny = len(lat) xmin, ymin, xmax, ymax = [long.min(), lat.min(), long.max(), lat.max()] xres = (xmax - xmin) / float(nx) yres = (ymax - ymin) / float(ny) geotransform = (xmin, xres, 0, ymax, 0, -yres)  #creates 1 raster band file dst_ds = gdal.getdriverbyname('gtiff').create('mygeotiff.tif', ny, nx, 1, gdal.gdt_float32)  dst_ds.setgeotransform(geotransform)    # specify coords srs = osr.spatialreference()            # establish encoding srs.importfromepsg(3857)                # wgs84 lat/long dst_ds.setprojection(srs.exporttowkt()) # export coords file dst_ds.getrasterband(1).writearray(x)   # write r-band raster dst_ds.flushcache()                     # write disk dst_ds = none                           # save, close

the process of geotiff creation above have largely sourced here:

how write/create geotiff rgb image file in python?

i have attempted based on understanding array want write raster 2d array 3 columns, 2 of co-ords , 1 of data. result check in snap black page white line along lhs.

so question follows:

how can extract necessary geotransform data netcdf file, adjust geotransform parameters appropriately , subsequently use extracted lat long + dataw arrays write geotiff file?


Comments

Post a Comment

Popular posts from this blog

commonjs - How to write a typescript definition file for a node module that exports a function? -

consider, toml node module can use: // toml.d.ts declare module toml { export function parse(value:string):any; } declare module "toml" { export = toml; } then: /// <reference path="../../../../../defs/toml/toml.d.ts"/> import toml = require('toml'); toml.parse(...); however, node module export single function, such 'glob' ( https://github.com/isaacs/node-glob ). the node usage of module be: var glob = require("glob") glob("*.js", {}, callback(err, files) { ... }); you'd naively expect this: // glob.d.ts declare function globs(paths:string, options:any, callback:{(err:any, files:string[]):void; ...but because typescripts 'import' semantics bit strange, seems can use 'import .. = require()' statement alias modules . trying call: /// <reference path="../../../../../defs/glob/glob.d.ts"/> import blog = require('glob'); results in: error ts2072: mod...
Read more

openid - Okta: Failed to get authorization code through API call -

i'm integrating okta own idp server using okta's api. i'm implementing authorization code flow following steps below: in own server, use /api/v1/authn endpoint sessiontoken. use sessiontoken obtain authorization calling endpoint: /oauth2/v1/authorize?client_id=" + clientid + "&sessiontoken=" + sessiontoken + "&response_type=code&response_mode=query&scope=openid&redirect_uri=" + redirecturl + "&state=evanyang&nonce=" it's supposed return response status code 302 , location header containing redirect url code value. however, keep getting response status code 200 , without location header, html body saying "you using unsupported browser." , "javascript disabled on browser." according api documentation: http://developer.okta.com/docs/api/resources/oidc.html#authentication-request , sessiontoken parameter sufficient this: an okta one-time sessiontoken. allows api-based ...
Read more

ios - Change Storyboard View using Seague -

i have following code checks connection server, , depending on want change views if no connection found. nsstring *connect = [nsstring stringwithcontentsofurl:[nsurl urlwithstring:@"http://myserver.com"] encoding:nsutf8stringencoding error:nil]; if (connect == null) { [self performseguewithidentifier:@"network_failure" sender:self]; } else { // other code here } my problem not change views. have added storyboard segue between 2 views identifier of network_failure . not have navigation controller or that, problem? thanks help! no not problem. first of all, in objectivec exists nil object oriented null . must use that. so: if(connect == nil) or better: if(!connect) so due fact nil not same of null , probably code does't enter in if scope . you have been able discover breakpoint. anyway: be sure set trigger segue (line) 2 view controller , , name network_failure (that not best name. should instead view1toview2 ) as, m...
Read more