Changeset 9741
- Timestamp:
- Oct 24, 2006, 2:29:31 PM (20 years ago)
- Location:
- trunk/icd-demo/src
- Files:
-
- 2 edited
-
icd-demo.c (modified) (3 diffs)
-
icd-test.c (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/icd-demo/src/icd-demo.c
r9127 r9741 11 11 12 12 if (argc != 3) { 13 14 13 fprintf (stderr, "USAGE: icd-demo (input) (output)\n"); 14 exit (2); 15 15 } 16 16 … … 38 38 while ((item = psListGetAndIncrement (iter)) != NULL) { 39 39 char *file = item->data.V; 40 40 parseFile (fits, file); 41 41 } 42 42 psFitsClose (fits); … … 58 58 psMetadata *blank = psMetadataLookupPtr (NULL, fd, "BLANK"); 59 59 if (blank) { 60 61 60 fprintf (stderr, "adding BLANK\n"); 61 psFitsWriteBlank (fits, blank, ""); 62 62 } 63 63 64 64 psMetadata *fitstable = psMetadataLookupPtr (NULL, fd, "FITSTABLE"); 65 65 if (fitstable) { 66 67 68 69 70 71 72 psArray *table = psArrayAlloc(10);73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 66 fprintf (stderr, "adding TABLE\n"); 67 psMetadata *header = psMetadataLookupPtr (NULL, fitstable, "HEADER"); 68 psMetadata *columns = psMetadataLookupPtr (NULL, fitstable, "COLUMNS"); 69 psMetadataIterator *iter = psMetadataIteratorAlloc (columns, PS_LIST_HEAD, NULL); 70 char *extname = psMetadataLookupStr (NULL, header, "EXTNAME"); 71 72 psArray *table = psArrayAllocEmpty (10); 73 74 // XXX insert unit values: we need a method to do this in psFitsWriteTable 75 int nCol = 1; 76 while ((col = psMetadataGetAndIncrement (iter)) != NULL) { 77 psMetadata *md = col->data.V; 78 char *unit = psMetadataLookupStr (NULL, md, "TUNIT"); 79 if (strcasecmp (unit, "NONE")) { 80 sprintf (unitWord, "TUNIT%d", nCol); 81 psMetadataAddStr (header, PS_LIST_TAIL, unitWord, 0, "", unit); 82 } 83 nCol ++; 84 } 85 psMetadataIteratorSet (iter, PS_LIST_HEAD); 86 87 // make a total of NN fake data columns 88 for (int i = 0; i < 10; i++) { 89 90 int nCol = 1; 91 92 psMetadata *row = psMetadataAlloc (); 93 94 while ((col = psMetadataGetAndIncrement (iter)) != NULL) { 95 96 psMetadata *md = col->data.V; 97 98 char *format = psMetadataLookupStr (NULL, md, "TFORM"); 99 char *comment = psMetadataLookupStr (NULL, md, "COMMENT"); 100 char *unit = psMetadataLookupStr (NULL, md, "TUNIT"); 101 if (strcasecmp (unit, "NONE")) { 102 sprintf (unitWord, "TUNIT%d", nCol); 103 psMetadataAddStr (header, PS_LIST_TAIL, unitWord, 0, "", unit); 104 } 105 106 if (!strcasecmp (format, "I")) { 107 int value; 108 109 char *minWord = psMetadataLookupStr (NULL, md, "TMIN"); 110 char *maxWord = psMetadataLookupStr (NULL, md, "TMAX"); 111 112 min = atoi (minWord); 113 if (!strcasecmp (minWord, "NONE")) min = 0; 114 max = atoi (maxWord); 115 if (!strcasecmp (maxWord, "NONE")) max = 100; 116 range = max - min; 117 118 rvalue = psRandomUniform (rnd); 119 value = range * rvalue + min; 120 121 psMetadataAdd (row, PS_LIST_TAIL, col->name, PS_DATA_S16, comment, value); 122 } 123 124 if (!strcasecmp (format, "J")) { 125 int value; 126 127 char *minWord = psMetadataLookupStr (NULL, md, "TMIN"); 128 char *maxWord = psMetadataLookupStr (NULL, md, "TMAX"); 129 130 min = atoi (minWord); 131 if (!strcasecmp (minWord, "NONE")) min = 0; 132 max = atoi (maxWord); 133 if (!strcasecmp (maxWord, "NONE")) max = 100; 134 range = max - min; 135 136 rvalue = psRandomUniform (rnd); 137 value = range * rvalue + min; 138 139 psMetadataAdd (row, PS_LIST_TAIL, col->name, PS_DATA_S32, comment, value); 140 } 141 142 if (!strcasecmp (format, "K")) { 143 long int value; 144 145 char *minWord = psMetadataLookupStr (NULL, md, "TMIN"); 146 char *maxWord = psMetadataLookupStr (NULL, md, "TMAX"); 147 148 min = atoi (minWord); 149 if (!strcasecmp (minWord, "NONE")) min = 0; 150 max = atoi (maxWord); 151 if (!strcasecmp (maxWord, "NONE")) max = 100; 152 range = max - min; 153 154 rvalue = psRandomUniform (rnd); 155 value = range * rvalue + min; 156 157 psMetadataAdd (row, PS_LIST_TAIL, col->name, PS_DATA_S64, comment, value); 158 } 159 160 if (!strcasecmp (format, "E")) { 161 float value; 162 163 char *minWord = psMetadataLookupStr (NULL, md, "TMIN"); 164 char *maxWord = psMetadataLookupStr (NULL, md, "TMAX"); 165 166 min = atoi (minWord); 167 if (!strcasecmp (minWord, "NONE")) min = 0; 168 max = atoi (maxWord); 169 if (!strcasecmp (maxWord, "NONE")) max = 100; 170 range = max - min; 171 172 rvalue = psRandomUniform (rnd); 173 value = range * rvalue + min; 174 175 psMetadataAdd (row, PS_LIST_TAIL, col->name, PS_DATA_F32, comment, value); 176 } 177 178 if (!strcasecmp (format, "D")) { 179 double value; 180 181 char *minWord = psMetadataLookupStr (NULL, md, "TMIN"); 182 char *maxWord = psMetadataLookupStr (NULL, md, "TMAX"); 183 184 min = atoi (minWord); 185 if (!strcasecmp (minWord, "NONE")) min = 0; 186 max = atoi (maxWord); 187 if (!strcasecmp (maxWord, "NONE")) max = 100; 188 range = max - min; 189 190 rvalue = psRandomUniform (rnd); 191 value = range * rvalue + min; 192 193 psMetadataAdd (row, PS_LIST_TAIL, col->name, PS_DATA_F64, comment, value); 194 } 195 196 if (!strcasecmp (format, "A")) { 197 char *value = NULL; 198 psStringAppend (&value, "word %03d", i); 199 psMetadataAdd (row, PS_LIST_TAIL, col->name, PS_DATA_STRING, "", value); 200 } 201 } 202 psArrayAdd (table, 100, row); 203 psFree (row); 204 psMetadataIteratorSet (iter, PS_LIST_HEAD); 205 } 206 psFitsInsertTable (fits, header, table, extname, true); 207 207 } 208 208 return 1; -
trunk/icd-demo/src/icd-test.c
r9127 r9741 7 7 8 8 if (argc != 3) { 9 10 9 fprintf (stderr, "USAGE: icd-demo (input) (output)\n"); 10 exit (2); 11 11 } 12 12 … … 15 15 16 16 { 17 18 19 17 psMetadata *header = psMetadataAlloc (); 18 psMetadataAddStr (header, PS_LIST_HEAD, "EXTID", 0, "", "blank"); 19 psFitsWriteBlank (fits, header, "blank"); 20 20 } 21 21 22 22 { 23 24 25 23 psMetadata *header = psMetadataAlloc (); 24 psMetadataAddStr (header, PS_LIST_HEAD, "EXTID", 0, "", "blank1"); 25 psFitsWriteBlank (fits, header, "blank1"); 26 26 } 27 27 28 28 { 29 30 29 psMetadata *header = psMetadataAlloc (); 30 psMetadataAddStr (header, PS_LIST_HEAD, "EXTID", 0, "", "table1"); 31 31 32 psArray *table = psArrayAlloc (100); 33 for (int i = 0; i < 10; i++) { 34 35 psMetadata *row = psMetadataAlloc (); 36 psMetadataAdd (row, PS_LIST_TAIL, "N", PS_DATA_S16, "", i); 37 psMetadataAdd (row, PS_LIST_TAIL, "F", PS_DATA_F32, "", i*0.2); 38 psMetadataAdd (row, PS_LIST_TAIL, "G", PS_DATA_F32, "", i*0.5); 39 psArrayAdd (table, 100, row); 40 psFree (row); 41 } 32 psArray *table = psArrayAllocEmpty (100); 33 for (int i = 0; i < 10; i++) { 42 34 43 psFitsInsertTable (fits, header, table, "table1", true); 35 psMetadata *row = psMetadataAlloc (); 36 psMetadataAdd (row, PS_LIST_TAIL, "N", PS_DATA_S16, "", i); 37 psMetadataAdd (row, PS_LIST_TAIL, "F", PS_DATA_F32, "", i*0.2); 38 psMetadataAdd (row, PS_LIST_TAIL, "G", PS_DATA_F32, "", i*0.5); 39 psArrayAdd (table, 100, row); 40 psFree (row); 41 } 42 43 psFitsInsertTable (fits, header, table, "table1", true); 44 44 } 45 45 46 46 { 47 48 49 47 psMetadata *header = psMetadataAlloc (); 48 psMetadataAddStr (header, PS_LIST_HEAD, "EXTID", 0, "", "blank2"); 49 psFitsWriteBlank (fits, header, "blank2"); 50 50 } 51 51 52 52 { 53 54 53 psMetadata *header = psMetadataAlloc (); 54 psMetadataAddStr (header, PS_LIST_HEAD, "EXTID", 0, "", "table2"); 55 55 56 psArray *table = psArrayAlloc (100); 57 for (int i = 0; i < 10; i++) { 58 59 psMetadata *row = psMetadataAlloc (); 60 psMetadataAdd (row, PS_LIST_TAIL, "N", PS_DATA_S16, "", i); 61 psMetadataAdd (row, PS_LIST_TAIL, "F", PS_DATA_F32, "", i*0.2); 62 psMetadataAdd (row, PS_LIST_TAIL, "G", PS_DATA_F32, "", i*0.5); 63 psArrayAdd (table, 100, row); 64 psFree (row); 65 } 56 psArray *table = psArrayAllocEmpty (100); 57 for (int i = 0; i < 10; i++) { 66 58 67 psFitsInsertTable (fits, header, table, "table2", true); 59 psMetadata *row = psMetadataAlloc (); 60 psMetadataAdd (row, PS_LIST_TAIL, "N", PS_DATA_S16, "", i); 61 psMetadataAdd (row, PS_LIST_TAIL, "F", PS_DATA_F32, "", i*0.2); 62 psMetadataAdd (row, PS_LIST_TAIL, "G", PS_DATA_F32, "", i*0.5); 63 psArrayAdd (table, 100, row); 64 psFree (row); 65 } 66 67 psFitsInsertTable (fits, header, table, "table2", true); 68 68 } 69 69
Note:
See TracChangeset
for help on using the changeset viewer.
