Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1228 1776 (N56W07 L= 87)  1041    233        1 (1130/1743-1334/1809)  -2.6 E04-W19 N54-N58
      2  625 1573 (N30E37 L= 43)  2045    403        3 ( 522/1471- 731/1641)  23.2 E44-E28 N23-N35
      3 1419 1537 (N38W23 L=104)   588     98        1 (1374/1508-1472/1558)  24.0 W19-W28 N35-N40
      4  386 1495 (N22E55 L= 26)   946    108        4 ( 360/1449- 414/1549) -60.3 E59-E51 N18-N25
      5 1791 1355 (N30W58 L=139)  1895    208        3 (1688/1324-1887/1400)  -1.8 W46-W77 N27-N32
      6  303 1205 (N03E52 L= 28)   711    134        2 ( 247/1149- 348/1244)  -5.2 E58-E47 N00-N05
      7  695  968 (S06E20 L= 61)   658    181        0 ( 638/ 899- 750/1046)  60.9 E22-E16 S10-S01
      8 1308  962 (N00W18 L= 98)   817    181        2 (1246/ 887-1348/1033)  63.4 W14-W20 S04-N04
      9 1821  659 (S13W68 L=148)   865    134        0 (1795/ 606-1853/ 716)  76.0 W65-W70 S16-S08
     10 1525  633 (S18W39 L=119)   536    217        1 (1448/ 596-1580/ 681)  24.0 W32-W42 S20-S14
     11 1577  521 (S24W48 L=128)   523    112        0 (1529/ 503-1633/ 553)  29.7 W43-W51 S25-S21
     12 1276  343 (S41W31 L=112)  1356    240        4 (1146/ 318-1372/ 382)   6.0 W20-W39 S43-S37
  

    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