Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1164 1669 (N45W25 L=235)   808    163        1 (1105/1629-1244/1699) -45.0 W21-W29 N40-N48
      2  636 1584 (N48E23 L=186)  1607    390        3 ( 481/1495- 788/1641) -18.4 E42-E09 N39-N53
      3 1605 1537 (N27W54 L=264)   752    122        3 (1542/1504-1634/1573) -29.7 W48-W59 N24-N28
      4 1120 1504 (N34W15 L=225)  1354    296        5 ( 988/1477-1273/1529) -11.3 W05-W25 N30-N37
      5 1764 1426 (N16W64 L=274)  1408    130        1 (1711/1375-1803/1461) -54.5 W59-W69 N12-N19
      6 1778 1267 (N07W58 L=267)   336     33        0 (1754/1247-1795/1285)   0.0 W55-W60 N05-N08
      7  511 1016 (N12E33 L=176)   351    107        1 ( 464/ 987- 564/1035)   8.1 E36-E29 N10-N12
      8  330  967 (N11E48 L=161)   282     75        0 ( 311/ 942- 358/1004) -63.4 E49-E45 N08-N12
      9  287  815 (N01E55 L=155)  1900    487        6 ( 122/ 765- 480/ 905) -11.9 E77-E38 S02-N07
     10  791  866 (S00E16 L=193)   331     85        1 ( 757/ 843- 826/ 901) -14.0 E18-E14 S01-N02
     11 1293  683 (S19W12 L=221)  2246    487        7 (1207/ 575-1416/ 881)  63.4 W05-W19 S26-S06
     12 1136  822 (S07W04 L=213)   341     82        2 (1106/ 790-1160/ 866)  36.9 W01-W05 S09-S05
     13 1096  739 (S12W00 L=210)   229     79        0 (1064/ 721-1134/ 761)   0.0 E01-W02 S13-S10
     14 1166  700 (S16W04 L=214)   366    180        0 (1090/ 678-1223/ 752)   7.1 W00-W08 S16-S12
     15 1556  717 (S21W31 L=240)   463     88        1 (1514/ 690-1597/ 742)  18.4 W27-W34 S22-S20
     16  571  645 (S11E35 L=175)   288     74        1 ( 550/ 618- 587/ 678)  45.0 E36-E33 S12-S09
     17  477  534 (S18E46 L=163)  1033    262        3 ( 373/ 469- 528/ 633) -63.4 E52-E42 S22-S10
     18  713  401 (S30E33 L=176)   531    129        0 ( 676/ 350- 745/ 454)  35.0 E38-E28 S33-S26
  

    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