Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1139 1799 (N57E04 L=101)   743    183        1 (1057/1774-1236/1818)   0.0 E13-W06 N54-N57
      2  430 1619 (N30E59 L= 46)  1611    295        1 ( 377/1498- 479/1694) -47.1 E71-E52 N21-N35
      3 1073 1614 (N39E05 L=100)  1085    323        2 ( 913/1576-1206/1676)  -2.4 E19-W05 N36-N41
      4  343 1589 (N27E68 L= 37)   816    100        1 ( 318/1543- 371/1631) -56.3 E71-E63 N23-N29
      5 1638 1368 (N30W41 L=146)  2494    445        5 (1437/1299-1821/1440)  -4.6 W23-W60 N26-N32
      6 1812 1108 (N15W56 L=161)   456    101        0 (1795/1064-1847/1134)  63.4 W54-W60 N12-N16
      7  368 1073 (S05E44 L= 61)   350     95        0 ( 337/1022- 400/1112)  63.4 E46-E41 S07-S01
      8  904 1030 (S01E07 L= 98)   408    113        1 ( 864/ 980- 944/1071)  68.2 E09-E04 S03-N02
      9  614 1003 (S06E25 L= 80)   254     98        0 ( 596/ 951- 634/1038)  78.7 E26-E23 S08-S03
     10  333  900 (S16E46 L= 59)   300     88        0 ( 315/ 852- 345/ 944) -71.6 E47-E44 S18-S12
     11 1188  708 (S17W14 L=119)   646    266        1 (1110/ 660-1244/ 772)  36.9 W09-W17 S20-S12
     12 1250  581 (S24W21 L=126)   281    101        0 (1206/ 564-1295/ 604)  38.7 W18-W23 S25-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