Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1301 1457 (N30W22 L=299)   448    117        0 (1283/1403-1319/1512)  76.0 W20-W24 N26-N34
      2  461 1381 (N28E41 L=236)   456    158        2 ( 391/1351- 530/1402)  -8.1 E48-E34 N26-N29
      3 1388 1360 (N23W27 L=304)   467     65        0 (1351/1337-1433/1381) -31.0 W24-W29 N21-N25
      4  460 1210 (N17E38 L=239)   266    100        1 ( 436/1156- 486/1249)  78.7 E39-E36 N14-N19
      5  560 1155 (N14E30 L=247)   238    113        0 ( 523/1119- 606/1177)  26.6 E32-E26 N11-N15
      6 1190  893 (S04W10 L=287)   238     45        1 (1163/ 867-1219/ 911)  18.4 W08-W11 S05-S03
      7 1742  758 (S17W52 L=329)  2548    407        5 (1661/ 611-1803/ 885) -69.8 W44-W59 S27-S08
      8 1295  832 (S09W16 L=293)   462     88        3 (1265/ 791-1320/ 877)  71.6 W14-W17 S11-S05
      9 1625  799 (S13W40 L=317)   484    179        1 (1587/ 726-1646/ 875)  71.6 W37-W41 S17-S08
     10 1823  683 (S23W65 L=342)  2479    217        4 (1778/ 586-1865/ 762)  79.7 W61-W70 S30-S18
     11  304  610 (S21E59 L=218)   783    137        2 ( 287/ 551- 319/ 680)  59.0 E62-E54 S25-S15
  

    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