Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  626 1727 (N64E22 L= 76)   383     67        2 ( 594/1710- 669/1743)  -9.5 E28-E15 N62-N66
      2 1716 1457 (N15W61 L=159)   431     58        2 (1691/1433-1741/1484) -56.3 W58-W63 N13-N17
      3  794 1415 (N35E07 L= 91)   445    162        0 ( 728/1362- 852/1440) -41.2 E11-E03 N30-N38
      4 1902 1116 (S09W69 L=166)   632     46        1 (1884/1090-1918/1141)  18.4 W65-W72 S11-S08
      5  771 1100 (N16E14 L= 84)   262    110        0 ( 732/1070- 819/1136) -51.3 E15-E11 N13-N18
      6 1812  967 (S15W53 L=151)   378     77        0 (1792/ 926-1828/1001)  36.9 W50-W55 S16-S13
      7 1330  900 (S07W16 L=113)   278     50        0 (1313/ 876-1357/ 930) -76.0 W14-W17 S08-S04
      8 1441  730 (S19W20 L=117)   632    242        0 (1382/ 647-1464/ 838)  51.3 W14-W22 S23-S13
      9 1752  550 (S43W51 L=149)   547     51        0 (1720/ 536-1783/ 566)  -5.7 W46-W56 S44-S41
     10  611  461 (S20E40 L= 57)   320     77        0 ( 573/ 448- 650/ 480)   0.0 E43-E36 S21-S19
     11  702  261 (S37E47 L= 51)  1851    225        1 ( 665/ 160- 753/ 354)  25.3 E71-E35 S47-S29
  

    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