Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1519 1636 (N35W41 L=126)   326     60        0 (1498/1614-1540/1659)  45.0 W38-W42 N33-N37
      2 1301 1508 (N27W20 L=106)   339     43        0 (1268/1480-1334/1533) -36.9 W18-W22 N24-N28
      3  703 1454 (N25E20 L= 66)   473    151        0 ( 638/1435- 792/1469)   0.0 E24-E13 N23-N26
      4 1103 1372 (N18W06 L= 92)   810    175        1 (1009/1322-1192/1402) -24.4 W00-W11 N14-N20
      5  561 1354 (N19E29 L= 56)   913    215        0 ( 470/1330- 666/1383)   0.0 E35-E21 N17-N21
      6 1459 1315 (N14W29 L=114)   472    149        0 (1408/1268-1525/1341)  12.5 W24-W33 N10-N15
      7  483 1142 (N07E34 L= 52)  2420    490        6 ( 408/ 980- 641/1331) -54.5 E40-E23 S03-N18
      8  288  652 (S22E56 L= 30)   770    132        1 ( 244/ 593- 332/ 702) -45.0 E58-E52 S25-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