OOPS (Re: missing points break up surface)

• To: visad-list@xxxxxxxxxxxxx
• Subject: OOPS (Re: missing points break up surface)
• From: "William L. Hibbard" <WHIBBARD@xxxxxxxxxxxxx>
• Date: Fri, 08 Jan 1999 07:26 CST

There is an error in my code for reconstructing an irregular
topology using the non-missing points.  A more accurate version
is:
 
  float[][] set_samples = domainSet.getSamples();
  // count non-missing samples
  int len = 0;
  for (int i=0; i<LengthX*LengthY; i++) {
    if (samples[0][i] == samples[0][i]) len++; // not missing
  }
  FlatField new_ff = null;
  if (len > 0) {
    float[][] new_set_samples = new float[2][len];
    double[][] new_samples = new double[2][len];
    len = 0;
    for (int i=0; i<LengthX*LengthY; i++) {
      if (samples[0][i] == samples[0][i]) {
        new_set_samples[0][len] = set_samples[0][i];
        new_set_samples[1][len] = set_samples[1][i];
        new_samples[0][len] = samples[0][i];
        new_samples[1][len] = samples[1][i];
        len++;
      }
    }
    Set new_domainSet = new Irregular2DSet( Domain, new_set_samples );
    new_ff = new FlatField( overfunc, new_domainSet );
    new_ff.setSamples(new_samples);
  }
  // now new_ff can be used in place of ff
  // if all points are missing, new_ff == null
 
> Since I don't have your data or ArcGridReader I was not able to test
> my code examples, so there may be a few bugs.
 
I was right.
 
----------------------------------------------------------
Bill Hibbard, SSEC, 1225 W. Dayton St., Madison, WI  53706
whibbard@xxxxxxxxxxxxx  608-263-4427  fax: 608-263-6738
http://www.ssec.wisc.edu/~billh/vis.html