Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1152 1578 (N35W26 L=157)  1840    411        6 ( 979/1522-1354/1613) -15.9 W13-W41 N28-N38
      2  913 1451 (N34W05 L=136)   471     96        0 ( 888/1409- 936/1503)  39.8 W02-W08 N31-N37
      3 1646  997 (S12W37 L=167)   312     73        0 (1613/ 969-1679/1030) -59.0 W35-W39 S14-S08
      4 1908  992 (S22W65 L=195)   609      7        0 (1880/ 961-1932/1030)  12.5 W60-W70 S22-S20
      5 1833  931 (S22W53 L=184)  1123    124        0 (1788/ 901-1882/ 963) -26.6 W48-W60 S24-S18
      6 1572  674 (S29W26 L=156)   685    116        2 (1509/ 644-1629/ 706) -16.7 W20-W30 S30-S26
      7 1816  657 (S40W53 L=183)   722    100        0 (1782/ 621-1848/ 702)  20.6 W48-W58 S42-S38
      8 1482  562 (S33W17 L=147)   801    136        1 (1436/ 514-1517/ 627)  32.0 W12-W20 S36-S29
      9  488  420 (S17E55 L= 76)   533     98        1 ( 434/ 389- 530/ 450) -56.3 E58-E51 S19-S13
     10  375  387 (S17E72 L= 58)  1562     75        1 ( 342/ 352- 403/ 430)  50.2 E79-E67 S20-S14
  

    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