Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1035 1755 (N37E23 L=100)  1132    189        2 ( 973/1693-1120/1790)  70.0 E26-E19 N30-N41
      2 1150 1567 (N28E08 L=116)   233     98        0 (1137/1520-1169/1615) -81.9 E09-E06 N24-N31
      3 1498 1497 (N34W16 L=140)   512     81        0 (1471/1459-1524/1537) -45.0 W13-W18 N31-N36
      4 1680 1302 (N28W34 L=158)   593    132        0 (1644/1247-1712/1359) -35.5 W30-W37 N25-N30
      5  960 1326 (N08E11 L=112)   587    181        2 ( 866/1296-1045/1350)  26.6 E16-E06 N05-N10
      6  489 1221 (S09E37 L= 87)   407    130        0 ( 452/1168- 514/1283)  90.0 E38-E35 S12-S05
      7  376 1178 (S13E44 L= 79)   272     52        0 ( 362/1143- 389/1208) -63.4 E45-E42 S15-S11
      8  905 1099 (S06E08 L=115)   603    137        2 ( 860/1039- 944/1158) -35.5 E11-E04 S08-S03
      9  639 1105 (S12E24 L=100)   410     84        1 ( 608/1071- 675/1140)  59.0 E25-E22 S14-S09
     10  818 1052 (S11E12 L=112)   228     36        0 ( 791/1030- 841/1073)  56.3 E13-E10 S12-S09
     11  885  989 (S12E07 L=117)   187     71        2 ( 855/ 968- 917/1025)  90.0 E08-E04 S14-S10
     12 1816  952 (N12W53 L=176)   601    128        1 (1769/ 918-1867/ 978)  36.9 W49-W58 N08-N14
     13  954  899 (S16E01 L=123)   444    153        2 ( 896/ 866-1031/ 937)  45.0 E03-W03 S18-S12
     14 1810  821 (N05W57 L=180)  1118    185        3 (1767/ 751-1853/ 896) -15.9 W49-W63 N01-N08
     15 1654  600 (S11W53 L=177)   507    100        2 (1625/ 550-1688/ 643)  90.0 W50-W54 S14-S07
     16 1122  559 (S30W20 L=144)   230     47        1 (1087/ 541-1143/ 579)  14.0 W17-W21 S30-S28
     17 1654  480 (S17W61 L=185)  1263    163        3 (1561/ 448-1720/ 506)  28.3 W54-W68 S21-S13
     18 1415  377 (S30W50 L=173)   587     56        1 (1380/ 349-1450/ 408)  53.1 W47-W52 S32-S28
  

    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