Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  176 1125 (S13E63 L=259)   847    146        0 ( 138/1086- 221/1152)  26.6 E66-E57 S15-S10
      2 1774 1063 (N12W48 L= 10)   779    241        0 (1703/ 953-1802/1137) -42.0 W41-W51 N05-N14
      3 1849 1065 (N14W55 L= 18)   382     96        0 (1827/1025-1865/1103) -38.7 W52-W58 N11-N15
      4  345  980 (S20E44 L=278)  2138    670        0 ( 163/ 834- 573/1078)  25.2 E61-E27 S30-S13
      5 1895  851 (N03W67 L= 29)   472     73        0 (1877/ 819-1906/ 895) -76.0 W65-W69 N01-N05
      6  845  722 (S28E05 L=317)   202     38        0 ( 820/ 703- 870/ 742)  56.3 E06-E03 S29-S26
      7 1519  707 (S15W37 L=360)   259     46        0 (1492/ 685-1541/ 730)  63.4 W35-W38 S16-S12
      8  368  633 (S42E43 L=279)   885    228        3 ( 320/ 564- 420/ 704)  51.3 E48-E37 S46-S36
      9 1758  645 (S12W60 L= 22)   581    193        0 (1696/ 617-1799/ 660)   4.4 W52-W65 S13-S10
     10 1504  637 (S19W39 L=  1)   324    108        0 (1445/ 624-1549/ 656)  23.2 W34-W41 S20-S16
     11 1571  396 (S31W56 L= 18)  3389    536        5 (1434/ 356-1753/ 437)  26.6 W45-W71 S35-S22
  

    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