Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1306 1793 (N58E04 L=217)   815     89        1 (1282/1753-1334/1833) -73.3 E07-W00 N52-N62
      2 1052 1624 (N34E15 L=206)   749    243        1 (1002/1525-1122/1715) -33.7 E22-E07 N28-N39
      3 1563 1538 (N44W26 L=248)   671    114        1 (1526/1491-1597/1569) -38.7 W22-W29 N40-N46
      4 1827 1222 (N30W55 L=276)  1376    232        1 (1752/1130-1898/1298) -17.5 W45-W66 N26-N33
      5 1236  772 (S12W19 L=240)   231     71        1 (1215/ 738-1256/ 800) -45.0 W16-W19 S14-S10
      6  234  768 (S37E58 L=164)   731     72        3 ( 203/ 732- 265/ 797)  20.0 E63-E52 S39-S34
      7  291  594 (S48E55 L=166)  1101    106        2 ( 261/ 549- 335/ 637)  41.6 E61-E47 S52-S44
      8 1478  596 (S16W40 L=262)   243    112        0 (1427/ 574-1515/ 633)  56.3 W37-W41 S18-S12
      9 1357  428 (S29W41 L=262)  2102    403        4 (1164/ 369-1519/ 487)  38.7 W29-W49 S37-S21
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y