Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  421 1669 (N33E57 L=125)   692    124        3 ( 391/1616- 458/1712) -82.9 E60-E54 N28-N36
      2 1348 1483 (N26W17 L=199)   826    246        0 (1247/1396-1421/1540) -30.3 W10-W22 N20-N29
      3  509 1494 (N20E39 L=143)   295     49        1 ( 481/1475- 540/1516) -11.3 E41-E36 N19-N21
      4  869 1386 (N15E12 L=169)   713    207        3 ( 797/1338- 927/1446)  -6.3 E17-E08 N12-N19
      5  630 1316 (N10E27 L=155)  1087    336        3 ( 509/1192- 730/1421) -37.4 E36-E18 N02-N15
      6 1865 1081 (N07W59 L=241)  2761    430        4 (1769/ 951-1944/1240) -47.0 W49-W71 N00-N15
      7  253 1009 (S10E53 L=129)   318    126        0 ( 235/ 976- 280/1054) -76.0 E54-E50 S11-S07
      8  441  910 (S15E37 L=144)  1021    320        3 ( 389/ 795- 491/1023)  66.8 E41-E33 S22-S08
      9  386  935 (S14E41 L=140)   365     96        1 ( 366/ 904- 426/ 977) -45.0 E42-E38 S15-S11
     10  745  868 (S16E16 L=166)   938    367        1 ( 635/ 783- 826/ 920)  18.4 E23-E10 S21-S13
     11  966  859 (S15E02 L=179)   454    190        2 ( 888/ 818-1026/ 890)  29.7 E06-W01 S17-S13
     12 1658  772 (S14W44 L=226)   275     67        0 (1641/ 737-1672/ 813)  90.0 W43-W45 S15-S11
     13  443  684 (S29E39 L=143)   281     66        0 ( 426/ 654- 456/ 725)  59.0 E40-E37 S31-S26
     14 1847  626 (S19W71 L=253)   759     91        3 (1817/ 582-1880/ 663)  50.2 W68-W75 S22-S16
     15 1134  421 (S42W15 L=196)   441    135        0 (1093/ 385-1202/ 465) -21.8 W10-W21 S43-S38
     16 1600  437 (S35W54 L=236)   390     71        1 (1572/ 403-1621/ 464) -18.4 W49-W58 S36-S33
  

    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