csv - How to extract data from NetCDF file using C# and Microsoft Scientific Data Set? -


i trying extract data netcdf file using microsoft's scientific data set libraries (sds). i'm restricted using ms sds , c#. i'm not programmer trade i'm struggling basics working. begin i'd write simple script extract data , write .csv file. using introduction doc , codeplex tutorials. i've tried write simple c# console app reads file , writes out console or ideally .csv file.

using sds 1.3 command line can view contents of test file: 

d:\netcdf>sds list test.nc [2] acpr of type single (time:85) (south_north:213) (west_east:165) [1] times of type sbyte (time:85) (datestrlen:19) d:\netcdf>

my script looks this: 

using system; using system.io; using sds = microsoft.research.science.data; using microsoft.research.science.data.imperative;   namespace netcdfconsoleapp {     class program     {         static void main(string[] args)         {             /// gets path netcdf file used data source.             var dataset = sds.dataset.open("d:\\netcdf\\test.nc?openmode=readonly");              sbyte[,] times = dataset.getdata<sbyte[,]>("times");             //int32[,] times = dataset.getdata<int32[,]>("times");              //single[] longitudes = dataset.getdata<single[]>("west_east");             //var latitudes = dataset.getdata<single[]>("south_north");             single[,,] datavalues = dataset.getdata<single[,,]>("acpr");              (int itime = 50; itime < 60; itime++)             {                 (int ilongitude = 130; ilongitude < 150; ilongitude++)                  {                     (int ilatitude = 130; ilatitude < 140; ilatitude++)                     {                         // write output data                         float thisvalue = datavalues[itime,ilatitude,ilongitude];                          console.writeline(itime);                         console.writeline(ilatitude);                         console.writeline(ilongitude);                         console.writeline(thisvalue);                      }                  }             }         }     } }

if comment out var times... line runs. i'm struggling sds read time dimension. if use sbyte complains variable doesn't exist. if use int32 complains converting string. 

system.invalidoperationexception unhandled   hresult=-2146233079   message=requested variable not exist in data set   source=scientificdataset   stacktrace:        @ microsoft.research.science.data.imperative.datasetextensions.findvariable(dataset dataset, func`2 predicate)        @ microsoft.research.science.data.imperative.datasetextensions.getdata[d](dataset dataset, string variablename)        @ netcdfconsoleapp.program.main(string[] args) in \\gdc-fs01\user$\prm\visual studio 2015\projects\netcdfconsoleapp\program.cs:line 16        @ system.appdomain._nexecuteassembly(runtimeassembly assembly, string[] args)        @ system.appdomain.executeassembly(string assemblyfile, evidence assemblysecurity, string[] args)        @ microsoft.visualstudio.hostingprocess.hostproc.runusersassembly()        @ system.threading.executioncontext.runinternal(executioncontext executioncontext, contextcallback callback, object state, boolean preservesyncctx)        @ system.threading.executioncontext.run(executioncontext executioncontext, contextcallback callback, object state, boolean preservesyncctx)        @ system.threading.executioncontext.run(executioncontext executioncontext, contextcallback callback, object state)        @ system.threading.threadhelper.threadstart()   innerexception:

what missing?

i think have solved puzzle. problem 2nd variable empty. i'm not sure if accident or design.

i have found date of when data supposed start in meta data field. i've modified code retrieves this, , writes console. 

using system; using system.io; using sds = microsoft.research.science.data; using microsoft.research.science.data.imperative;   namespace netcdfconsoleapp {     class program     {         static void main(string[] args)         {             /// gets path netcdf file used data source.             var dataset = sds.dataset.open("d:\\netcdf\\test.nc?openmode=readonly");              string dt = (string)dataset.metadata["start_date"];              single[,,] datavalues = dataset.getdata<single[,,]>("acpr");              (int itime = 50; itime < 60; itime++)             {                 (int ilongitude = 130; ilongitude < 150; ilongitude++)                 {                     (int ilatitude = 130; ilatitude < 140; ilatitude++)                     {                         // write output data                         float thisvalue = datavalues[itime,ilatitude,ilongitude];                                                                         console.writeline(dt.tostring() + ',' + itime.tostring() + ',' + ilatitude.tostring() + ',' + ilongitude.tostring() + ',' + thisvalue.tostring());                      }                                  }             }                         console.readline();         }                } }

i've struggled i'm sharing in hope of use else.

one thing found useful discussion tab on codeplex has lots of useful code snippets.


Comments

Post a Comment

Popular posts from this blog

javascript - How to bind ViewModel Store to View? -

i'm pretty new ext js , trying embed multiselect inside panel . the viewmodel has stores property can see here: ext.define('test.view.controls.search.searchfiltermodel', { extend: 'ext.app.viewmodel', alias: 'viewmodel.filter', data: { title: '' }, stores: { test: { fields: [ 'id', 'name' ], proxy: { type: 'ajax', url: 'api/test', reader: 'array' }, autoload: true } } }); i bind in view this: viewmodel: { type: 'filter' }, layout: 'fit', border: 1, plain: true, scrollable: 'y', layout: 'fit', bind: { title: '{title}', }, items: { xtype: 'multiselect', scrollable: false, allowblank: true, ddreorder: true, bind: { store: '{test}' }, valuefield: 'id', displayfield: 'name' } ...
Read more

recursion - Can every recursive algorithm be improved with dynamic programming? -

i first year undergraduate csc student looking competitive programming. recursion involves defining , solving sub problems. understand, top down dynamic programming (dp) involves memoizing solutions sub problems reduce time complexity of algorithm. can top down dp used improve efficiency of every recursive algorithm overlapping sub problems? dp fail work , how can identify this? the short answer is: yes. however, there constraints. obvious 1 recursive calls must overlap. i.e. during execution of algorithm, recursive function must called multiple times same parameters. lets truncate recursion tree memoization. can use memoization reduce number of calls. however, reduction of calls comes price. need store results somewhere. next obvious constraint need have enough memory. comes not-so obvious constraint. memory access requires time. first need find result stored , maybe copy location. in cases, might faster let recursion calculate result instead of loading somewher...
Read more

c - Why does alarm() cause fgets() to stop waiting? -

i playing around signals in c. main function asks input using fgets(name, 30, stdin) , , sits there , waits. set alarm alarm(3) , , reassigned sigalrm call function myalarm calls system("say pay attention") . after alarm goes off, fgets() stops waiting input , main fn continues on. happens if change myalarm set variable , nothing it. void myalarm(int sig) { //system("say pay attention"); int x = 0; } int catch_signal(int sig, void (*handler)(int)) { // when signal comes in, "catch" , "handle it" in way want struct sigaction action; // create new sigaction action.sa_handler = handler; // set it's sa_handler attribute function specified in header sigemptyset(&action.sa_mask); // "turn signals in sa_mask off?" "set sa_mask contian no signals, i.e. nothing masked?" action.sa_flags = 0; // not sure, looks aren't using of available flags, whatever may return sigaction(sig, ...
Read more