Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  974 1765 (N38E26 L=271)   790    108        1 ( 934/1738-1035/1797)  38.7 E29-E21 N34-N39
      2  366 1611 (N14E63 L=234)  1893    331        3 ( 336/1543- 387/1689) -45.0 E69-E57 N10-N20
      3  455 1462 (N06E47 L=250)   805    198        3 ( 394/1417- 518/1511) -14.0 E53-E41 N04-N08
      4  882 1344 (N08E16 L=281)  1244    296        3 ( 851/1215- 908/1483) -61.9 E19-E11 N00-N16
      5 1688 1108 (N15W39 L=335)   935    152        1 (1649/1025-1720/1183) -45.0 W34-W42 N11-N18
      6  483 1088 (S15E34 L=262)   347    120        0 ( 429/1073- 539/1105)   7.1 E38-E30 S15-S14
      7  701 1070 (S12E20 L=277)   279     67        0 ( 663/1046- 740/1094)  -9.5 E22-E16 S12-S10
      8  227 1022 (S23E55 L=242)  1771    365        2 ( 143/ 949- 327/1085)  29.4 E62-E45 S28-S18
      9 1647  863 (N00W42 L=338)   303     77        0 (1613/ 845-1679/ 883)  26.6 W39-W44 S01-N01
     10 1726  546 (S13W62 L=359)  1153    130        0 (1666/ 516-1764/ 590)  45.0 W57-W66 S16-S10
  

    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