Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1617 1561 (N32W57 L=202)   651    112        0 (1571/1530-1671/1582) -16.7 W51-W61 N29-N33
      2 1148 1534 (N37W15 L=160)   343     79        0 (1101/1509-1187/1557)  15.9 W11-W18 N35-N38
      3 1569 1499 (N29W48 L=193)   688    115        1 (1539/1442-1597/1550)  82.9 W46-W49 N24-N32
      4 1732 1477 (N24W65 L=210)   913    113        1 (1697/1426-1775/1521) -90.0 W58-W68 N20-N27
      5  828 1357 (N28E10 L=135)   393    127        0 ( 800/1291- 870/1407) -63.4 E12-E07 N23-N31
      6  417 1287 (N26E42 L=103)  2449    419        6 ( 309/1174- 588/1341)  11.8 E52-E27 N18-N29
      7  160 1228 (N22E74 L= 71)  2490    162        3 ( 138/1150- 180/1295) -83.7 E84-E69 N17-N26
      8  416 1083 (N13E41 L=104)   538    100        1 ( 371/1051- 452/1112) -26.6 E44-E37 N10-N14
      9 1752  856 (S13W50 L=195)   615     45        0 (1722/ 827-1778/ 883)  11.3 W47-W53 S15-S11
     10  705  800 (S06E22 L=123)   244     92        0 ( 662/ 777- 742/ 816) -21.8 E24-E19 S07-S04
     11  266  741 (S08E59 L= 86)   539     70        0 ( 242/ 715- 303/ 760)  -9.5 E61-E55 S09-S06
     12 1083  479 (S31E01 L=144)   262     99        1 (1031/ 464-1129/ 494) -21.8 E04-W01 S32-S28
     13  460  402 (S32E55 L= 90)   669    112        0 ( 438/ 353- 491/ 461)  60.9 E57-E52 S36-S27
     14  884  405 (S34E17 L=128)   401    137        2 ( 823/ 389- 955/ 422)  -5.7 E21-E11 S35-S32
     15 1263  252 (S54W13 L=159)   425     37        0 (1241/ 236-1289/ 272) -26.6 W10-W16 S56-S51
     16  949  183 (S57E23 L=122)   795    106        4 ( 895/ 165-1007/ 203)   3.8 E30-E15 S59-S54
  

    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