Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1600 1684 (N24W67 L= 27)   482     69        0 (1557/1673-1634/1696) -20.6 W62-W71 N22-N25
      2 1414 1577 (N24W42 L=  2)   872    187        2 (1333/1539-1487/1602) -36.9 W37-W46 N20-N28
      3 1788 1472 (N07W69 L= 29)  3519    324       12 (1726/1398-1857/1581)  85.6 W60-W77 N01-N14
      4 1617 1486 (N13W52 L= 12)   319    148        2 (1569/1457-1650/1514) -26.6 W48-W55 N11-N15
      5 1284 1108 (N03W16 L=337)   219     64        0 (1267/1087-1309/1144)  71.6 W15-W17 N01-N04
      6  402 1094 (N23E38 L=282)   990    273        3 ( 304/1055- 511/1138)   0.0 E47-E28 N20-N25
      7 1626 1022 (S12W36 L=356)   698    205        6 (1583/ 947-1660/1079)  60.3 W32-W38 S15-S07
      8 1779  783 (S31W46 L=  6)  1553    188        2 (1748/ 713-1807/ 883)  68.2 W42-W49 S35-S24
      9  360  674 (N00E52 L=268)   299     80        2 ( 334/ 638- 394/ 691)   9.5 E54-E48 S01-N01
     10 1730  639 (S39W40 L=  0)   340     54        1 (1705/ 625-1759/ 662)  -8.1 W37-W44 S40-S36
     11  468  403 (S18E57 L=263)  2443    207        7 ( 395/ 331- 542/ 479) -31.0 E64-E51 S22-S12
     12  763  284 (S33E43 L=277)   413    107        1 ( 719/ 251- 809/ 313)  10.3 E49-E38 S34-S31
  

    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