Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1656 1460 (N20W54 L=185)   457    132        1 (1614/1434-1693/1476)   0.0 W48-W57 N19-N21
      2 1207 1201 (N14W14 L=146)   588    136        0 (1178/1131-1233/1272)  77.5 W12-W15 N09-N18
      3 1637 1232 (N08W44 L=176)   297     51        2 (1614/1209-1664/1257) -53.1 W42-W45 N05-N09
      4 1663 1039 (S04W42 L=173)   364     72        0 (1627/1012-1694/1067)   9.5 W38-W44 S05-S02
      5 1710  873 (S15W44 L=175)   766    131        3 (1694/ 804-1730/ 942)  77.5 W41-W45 S19-S10
      6 1922  898 (S20W68 L=200)   893    185        1 (1900/ 851-1942/ 942)  51.3 W64-W71 S22-S17
      7  347  653 (S08E54 L= 78)   902    160        4 ( 309/ 578- 393/ 741) -84.8 E56-E50 S13-S02
      8 1326  678 (S20W13 L=145)   359    111        1 (1312/ 622-1341/ 738)  76.0 W12-W15 S23-S15
      9 1605  637 (S28W34 L=166)   393     68        1 (1574/ 602-1633/ 667)  21.8 W31-W36 S30-S26
     10  822  514 (S22E22 L=110)  1187    265        2 ( 713/ 459- 943/ 573)  12.5 E30-E12 S25-S19
     11  455  506 (S19E49 L= 83)   360     47        0 ( 427/ 486- 483/ 526) -45.0 E51-E46 S20-S17
     12  589  327 (S33E47 L= 85)   833    155        2 ( 556/ 291- 619/ 361)  33.7 E51-E42 S36-S30
  

    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