Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1031 1876 (N63W42 L=134)  1106    185        2 ( 981/1823-1075/1905)  54.0 W36-W47 N57-N68
      2 1515 1450 (N19W43 L=135)   938    200        1 (1434/1422-1610/1486) -38.7 W38-W49 N15-N23
      3 1779 1248 (N01W57 L=148)  2233    305        3 (1684/1110-1828/1366) -86.4 W51-W61 S07-N09
      4 1110 1210 (N16W09 L=101)   409    134        0 (1055/1172-1143/1266)  23.2 W05-W12 N14-N18
      5  195 1101 (N26E59 L= 33)   590    128        0 ( 165/1073- 229/1135) -29.7 E62-E55 N24-N28
      6 1811  950 (S17W52 L=143)   908    136        2 (1760/ 914-1843/ 991) -18.4 W46-W55 S19-S15
      7  722  590 (S12E28 L= 64)  1274    352        0 ( 667/ 452- 806/ 722)  41.4 E36-E19 S20-S05
      8 1542  559 (S34W24 L=116)   523    121        1 (1497/ 530-1581/ 586) -45.0 W20-W27 S37-S31
      9 1568  501 (S40W27 L=118)  1055    173        3 (1502/ 459-1637/ 526) -21.0 W20-W33 S43-S36
     10  502  414 (S20E52 L= 40)   387    103        2 ( 476/ 388- 522/ 439)  31.0 E53-E48 S22-S18
  

    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