Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1003 1879 (N50E32 L=196)   585     54        2 ( 973/1858-1033/1904)  -6.3 E36-E27 N48-N53
      2  399 1452 (N06E51 L=177)   762    158        0 ( 355/1389- 452/1502)  69.4 E53-E48 N01-N09
      3  538 1465 (N09E40 L=188)   489    119        1 ( 482/1425- 596/1492)  -5.7 E45-E35 N07-N10
      4  284 1414 (N02E61 L=168)   531     79        2 ( 246/1388- 313/1439)  33.7 E64-E57 N00-N04
      5 1083 1400 (N16E05 L=223)   213     54        0 (1059/1374-1108/1426) -36.9 E06-E02 N14-N17
      6 1788 1335 (N32W48 L=276)   403    115        0 (1775/1283-1803/1394)  77.5 W46-W50 N27-N36
      7 1568 1203 (N16W29 L=257)   349    105        2 (1531/1175-1619/1234)  -8.1 W26-W33 N14-N17
      8  211  852 (S32E60 L=168)   367     40        0 ( 193/ 831- 229/ 874) -45.0 E62-E57 S33-S30
      9 1642  599 (S15W52 L=281)   631    157        1 (1573/ 538-1690/ 644)  60.9 W49-W55 S20-S11
     10 1487  360 (S33W56 L=285)   337     52        0 (1463/ 347-1517/ 381)   0.0 W52-W60 S33-S31
  

    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