Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1073 1828 (N48E17 L=195)  1335    297        3 (1033/1748-1117/1903) -74.5 E20-E12 N40-N58
      2  363 1509 (N13E55 L=157)   382     93        0 ( 330/1476- 402/1526)  45.0 E57-E51 N10-N14
      3 1446 1409 (N24W19 L=231)   685    116        1 (1412/1352-1484/1459)  66.0 W17-W21 N19-N28
      4  359 1378 (N04E51 L=161)   580    109        0 ( 337/1318- 388/1435)  82.9 E52-E49 N00-N08
      5 1740 1297 (N24W43 L=255)  1204    165        2 (1699/1216-1777/1386) -56.3 W39-W46 N20-N29
      6  794  967 (S14E13 L=199)  1133    322        6 ( 684/ 830- 882/1125)  68.2 E18-E09 S23-S03
      7 1624 1046 (N06W37 L=248)   426    113        1 (1577/1012-1665/1098) -23.2 W32-W39 N04-N08
      8  330  936 (S22E45 L=167)   993    324        1 ( 233/ 841- 411/1032)  42.9 E53-E39 S28-S15
      9  425  773 (S31E37 L=175)  1427    424        0 ( 298/ 688- 559/ 894)  37.3 E48-E26 S37-S21
     10 1626  656 (S15W47 L=259)  3994    935        4 (1247/ 589-1873/ 714)  18.8 W23-W76 S27-S07
  

    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