Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1734 1438 (N42W38 L= 77)   374     80        0 (1713/1418-1749/1467)  59.0 W35-W39 N39-N44
      2 1342 1355 (N21W09 L= 47)   309    127        0 (1293/1298-1385/1414) -32.0 W04-W12 N18-N23
      3 1150 1342 (N15E02 L= 37)   267     58        0 (1116/1322-1178/1367)  56.3 E03-W00 N13-N16
      4 1741 1107 (N20W42 L= 81)   973    169        0 (1718/1024-1762/1209) -59.0 W38-W45 N15-N25
      5 1189  967 (S05W11 L= 49)   311     61        1 (1151/ 947-1216/ 994)  45.0 W08-W12 S06-S03
      6 1870  926 (N15W60 L= 98)  1117    115        1 (1849/ 872-1892/ 986) -63.4 W56-W62 N11-N17
      7  522  824 (S31E26 L= 13)  1288    432        2 ( 336/ 771- 692/ 890)   5.9 E42-E13 S33-S27
      8 1350  766 (S11W26 L= 64)   534    143        1 (1316/ 698-1407/ 834)  90.0 W24-W27 S15-S05
      9 1238  618 (S22W25 L= 64)   436     76        1 (1191/ 589-1278/ 646)  51.3 W22-W27 S24-S19
     10 1612  351 (S23W66 L=105)  1227    142        1 (1553/ 324-1652/ 388)  54.5 W62-W69 S26-S19
  

    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