Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
smurf: Allow makemap to export compressed NOI model
Since skyloop is to support noi.box_size, it will need to cache the NOI model created by the first invocation of makemap. But this could be another fill size time-series data set since previously the NOI model could only be exprted as the Variance component of the RES model. Now, by setting noi.export, the NOI model can be exported in compressed form (i.e. only one value per noi.box_size time slices of raw data), without lots of unneeded extras like quality and residuals, etc. , but this
- Loading branch information
David Berry
committed
Sep 27, 2013
1 parent
a4e1407
commit dd0d74c
Showing
11 changed files
with
394 additions
and
121 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -0,0 +1,181 @@ | ||
| /* | ||
| *+ | ||
| * Name: | ||
| * smf_export_noi | ||
| * Purpose: | ||
| * Export the NOI model to an NDF. | ||
| * Language: | ||
| * Starlink ANSI C | ||
| * Type of Module: | ||
| * C function | ||
| * Invocation: | ||
| * void smf_export_noi( smfData *noi, const char *name, int boxsize, | ||
| * int *status ) | ||
| * Arguments: | ||
| * noi = smfData * (Given) | ||
| * The NOI model. | ||
| * name = const char * (Given) | ||
| * The name of the NDF to create. | ||
| * boxsize = int (Given) | ||
| * The number of samples in one noise box, corresponding to | ||
| * NOI.BOX_SIZE. | ||
| * status = int* (Given and Returned) | ||
| * Pointer to global status. | ||
| * Description: | ||
| * This function exports the compressed NOI model to a new NDF. It can | ||
| * be imported again via function smf_import_noi. | ||
| * Authors: | ||
| * David S Berry (JAC, Hawaii) | ||
| * {enter_new_authors_here} | ||
| * History: | ||
| * 24-SEP-2013 (DSB): | ||
| * Original version. | ||
| * {enter_further_changes_here} | ||
| * Copyright: | ||
| * Copyright (C) 2013 Science & Technology Facilities Council. | ||
| * All Rights Reserved. | ||
| * Licence: | ||
| * This program is free software; you can redistribute it and/or | ||
| * modify it under the terms of the GNU General Public License as | ||
| * published by the Free Software Foundation; either version 3 of | ||
| * the License, or (at your option) any later version. | ||
| * | ||
| * This program is distributed in the hope that it will be | ||
| * useful, but WITHOUT ANY WARRANTY; without even the implied | ||
| * warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR | ||
| * PURPOSE. See the GNU General Public License for more details. | ||
| * | ||
| * You should have received a copy of the GNU General Public | ||
| * License along with this program; if not, write to the Free | ||
| * Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, | ||
| * MA 02110-1301, USA | ||
| * Bugs: | ||
| * {note_any_bugs_here} | ||
| *- | ||
| */ | ||
|
|
||
| /* Starlink includes */ | ||
| #include "sae_par.h" | ||
| #include "mers.h" | ||
| #include "ndf.h" | ||
| #include "star/hds.h" | ||
| #include "ast.h" | ||
|
|
||
| /* SMURF includes */ | ||
| #include "libsmf/smf.h" | ||
|
|
||
| void smf_export_noi( smfData *noi, const char *name, int boxsize, int *status ){ | ||
|
|
||
| /* Local Variables */ | ||
| HDSLoc *xloc = NULL; | ||
| dim_t ntslice; | ||
| dim_t nrows; | ||
| dim_t ncols; | ||
| dim_t nbolo; | ||
| double *ip; | ||
| double *dataptr; | ||
| double *dp; | ||
| double *pd; | ||
| int el; | ||
| int ibolo; | ||
| int indf; | ||
| int itime; | ||
| int lbnd[ 3 ]; | ||
| int nz; | ||
| int place; | ||
| int ubnd[ 3 ]; | ||
| size_t bstride; | ||
| int iz; | ||
|
|
||
| /* Check inherited status. */ | ||
| if( *status != SAI__OK ) return; | ||
|
|
||
| /* Get the dimensions of the NOI model. */ | ||
| smf_get_dims( noi, &nrows, &ncols, &nbolo, &ntslice, NULL, &bstride, | ||
| NULL, status ); | ||
| if( ntslice == 1 ) boxsize = 0; | ||
|
|
||
| /* Determine the number of boxes to use. This is the length of 3rd axis. | ||
| of the new NDF. A boxsize of zero means "there is one box". */ | ||
| nz = 1; | ||
| if( boxsize > 0 ) { | ||
| nz = ntslice/boxsize; | ||
| if( nz == 0 ) nz = 1; | ||
| } else if( ntslice > 1 ){ | ||
| *status = SAI__ERROR; | ||
| errRepf("", "smf_export_noi: boxsize is zero but ntslice (%d) is " | ||
| "larger than 1 (programming error).", status, (int) ntslice ); | ||
| } | ||
|
|
||
| /* Get a pointer to the NOI data values. */ | ||
| dataptr = noi->pntr[ 0 ]; | ||
|
|
||
| /* Create the NDF. */ | ||
| ndfPlace( NULL, name, &place, status ); | ||
| lbnd[ 0 ] = 1; | ||
| lbnd[ 1 ] = 1; | ||
| lbnd[ 2 ] = 1; | ||
| ubnd[ 0 ] = ncols; | ||
| ubnd[ 1 ] = nrows; | ||
| ubnd[ 2 ] = nz; | ||
| ndfNew( "_DOUBLE", 3, lbnd, ubnd, &place, &indf, status ); | ||
|
|
||
| /* Map the Data array of the NDF and copy the NOI values into it. */ | ||
| ndfMap( indf, "DATA", "_DOUBLE", "WRITE", (void **) &ip, &el, status ); | ||
| if( *status == SAI__OK ) { | ||
|
|
||
| /* Initialise the time slice at the middle of the current box in the model. */ | ||
| itime = ( nz == 1 ) ? 0 : boxsize/2; | ||
|
|
||
| /* Loop round all planes in the NDF. */ | ||
| for( iz = 0; iz < nz; iz++ ) { | ||
|
|
||
| /* First deal with a time ordered model. The NOI model is filled with 1.0 | ||
| values by smf_model_create, and can also be set to zero to indicate | ||
| missing values. Therefore convert both these values into VAL__BADD. */ | ||
| if( bstride == 1 ){ | ||
| pd = ip; | ||
| dp = dataptr + itime*nbolo; | ||
| for( ibolo = 0; ibolo < (int) nbolo; ibolo++,dp++ ) { | ||
| *(pd++) = (*dp == 0.0 || *dp == 1.0) ? VAL__BADD : *dp; | ||
| } | ||
|
|
||
| /* Now deal with a bolo ordered model. The NOI model is filled with 1.0 | ||
| values by smf_model_create, and can also be set to zero to indicate | ||
| missing values. Therefore convert both these values into VAL__BADD. */ | ||
| } else { | ||
| pd = ip; | ||
| dp = dataptr + itime; | ||
| for( ibolo = 0; ibolo < (int) nbolo; ibolo++ ) { | ||
| *(pd++) = (*dp == 0.0 || *dp == 1.0) ? VAL__BADD : *dp; | ||
| dp += ntslice; | ||
| } | ||
| } | ||
|
|
||
| /* Point to the start of the next plane in the NDF. */ | ||
| ip += nbolo; | ||
| itime += boxsize; | ||
| } | ||
|
|
||
| /* Store the box size as an extension item in the NDF. */ | ||
| ndfXnew( indf, SMURF__EXTNAME, SMURF__EXTTYPE, 0, NULL, &xloc, status ); | ||
| ndfXpt0i( boxsize, indf, SMURF__EXTNAME, "NOI_BOXSIZE", status ); | ||
| datAnnul( &xloc, status ); | ||
| } | ||
|
|
||
| /* Annul the NDF identifier. */ | ||
| ndfAnnul( &indf, status ); | ||
|
|
||
| } | ||
|
|
Oops, something went wrong.