Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1347 1383 (N21W28 L=195)   666    225        0 (1248/1358-1451/1405) -21.0 W22-W35 N18-N24
      2 1574 1078 (S01W35 L=202)  1312    386        2 (1530/ 930-1640/1236)  86.8 W32-W40 S11-N07
      3 1698 1026 (S07W44 L=210)   329     74        0 (1685/ 992-1711/1064)  63.4 W42-W45 S09-S05
      4 1810  915 (S17W54 L=220)  1137    205        4 (1766/ 823-1841/ 992)  73.3 W48-W57 S22-S12
      5 1463  769 (S17W23 L=189)   553    207        0 (1425/ 689-1481/ 870)  59.0 W19-W25 S21-S11
      6 1065  821 (S06E01 L=165)   267    122        2 (1011/ 775-1106/ 857)  26.6 E05-W01 S07-S04
      7 1250  699 (S17W08 L=174)   521    210        3 (1165/ 658-1306/ 763)  24.0 W02-W11 S18-S13
      8 1678  592 (S35W41 L=207)  2595    432        3 (1615/ 472-1747/ 722)  58.0 W35-W49 S44-S27
      9  354  537 (S13E59 L=107)  1302    229        6 ( 330/ 450- 397/ 644)  69.0 E67-E54 S19-S06
     10 1160  408 (S35E04 L=163)  3817    598       17 ( 940/ 237-1325/ 599)  25.3 E25-W11 S46-S24
     11  569  545 (S15E40 L=127)   559    155        2 ( 500/ 514- 651/ 560)  -9.5 E44-E33 S16-S13
  

    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