Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  289 1572 (N12E74 L=277)  1579    120        1 ( 254/1532- 319/1617)  66.8 E78-E70 N08-N15
      2  223 1200 (S12E60 L=291)  1303    245        3 ( 176/1120- 278/1273)  73.3 E63-E55 S16-S06
      3 1822 1065 (N20W53 L= 43)  1181    249        2 (1759/1041-1896/1103)  11.3 W46-W61 N18-N22
      4 1409  993 (N04W23 L= 14)   424    105        3 (1385/ 938-1426/1050) -63.4 W21-W25 N01-N07
      5 1685  897 (N06W45 L= 35)   512    119        0 (1650/ 839-1716/ 959)  76.0 W43-W45 N02-N10
      6  883  676 (S27W01 L=352)   428    143        0 ( 794/ 638- 937/ 707)   0.0 E05-W05 S28-S26
      7 1724  577 (S10W62 L= 52)  1332    260        5 (1692/ 512-1762/ 631) -76.0 W56-W69 S13-S05
      8 1617  542 (S15W53 L= 44)   675     96        2 (1585/ 503-1656/ 587) -80.5 W50-W56 S17-S11
      9 1083  387 (S40W25 L= 16)  1788    390        2 ( 948/ 333-1247/ 493)   0.0 W11-W38 S43-S36
     10 1595  369 (S25W63 L= 53)  1166    258        0 (1535/ 352-1665/ 388)  15.9 W56-W70 S27-S21
  

    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