Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1330 1658 (N45W02 L=169)   378     81        0 (1286/1644-1367/1672)  26.6 E01-W05 N42-N46
      2  463 1599 (N16E56 L=111)  1217    206        5 ( 427/1516- 523/1664) -79.7 E61-E52 N10-N21
      3 1403 1549 (N39W10 L=177)   268     59        1 (1372/1528-1431/1569) -11.3 W07-W12 N37-N39
      4 1639 1306 (N30W33 L=200)  1915    380        2 (1494/1203-1780/1415) -10.5 W19-W48 N27-N32
      5 1736 1142 (N22W44 L=211)   398    141        1 (1692/1110-1776/1166)  35.5 W40-W48 N19-N24
      6 1051 1092 (N01E00 L=167)   387    112        1 (1014/1045-1079/1152) -38.7 E03-W02 S01-N03
      7 1624 1106 (N17W35 L=202)   319      0       -1 (1603/1075-1643/1132) -56.3 W33-W36 N15-N18
      8 1088  960 (S06W05 L=172)   489     97        3 (1052/ 921-1123/1005)  80.5 W03-W07 S08-S02
      9  314  900 (S28E46 L=121)   288     52        1 ( 302/ 868- 335/ 927) -76.0 E47-E44 S30-S26
     10  174  830 (S36E66 L=101)   625     98        2 ( 162/ 781- 187/ 884)  48.8 E70-E63 S39-S31
     11 1002  843 (S15W03 L=170)   333     77        0 ( 965/ 818-1041/ 865)  45.0 W01-W05 S16-S12
     12 1444  562 (S18W40 L=207)   274     36        1 (1424/ 543-1466/ 584)  90.0 W38-W41 S19-S16
  

    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