Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1273 1683 (N35W38 L=141)   648    173        0 (1187/1653-1354/1711) -45.0 W33-W42 N30-N39
      2 1695 1634 (N19W70 L=173)  1377    111        3 (1667/1581-1716/1687)  11.9 W64-W84 N16-N21
      3 1320 1567 (N26W34 L=137)   764    162        0 (1261/1500-1387/1606)   8.7 W27-W40 N23-N28
      4 1093 1408 (N23W14 L=117)   606    262        0 ( 989/1321-1169/1477) -70.3 W09-W16 N16-N30
      5  868 1284 (N22E03 L=100)   849    240        1 ( 812/1202- 956/1364) -76.0 E05-W00 N15-N27
      6 1752 1163 (S07W49 L=152)  2359    400        1 (1696/ 972-1794/1317)  59.5 W41-W53 S16-N02
      7  613 1219 (N24E20 L= 83)  2285    582        5 ( 369/1073- 806/1297) -32.7 E38-E09 N11-N32
      8  493 1086 (N19E31 L= 72)   271    113        0 ( 438/1071- 547/1102) -29.7 E34-E27 N17-N21
      9 1850  866 (S27W54 L=157)  2336    229        4 (1741/ 807-1932/ 934)   5.0 W41-W67 S29-S24
     10 1557  846 (S19W27 L=130)  1560    332        3 (1428/ 736-1685/ 920)  10.3 W16-W38 S22-S16
     11  394  583 (S07E52 L= 51)  4748    457       16 ( 303/ 447- 468/ 828)  52.7 E63-E42 S15-N06
     12  765  722 (S07E22 L= 81)   408    118        0 ( 728/ 677- 788/ 774)  45.0 E25-E19 S08-S03
     13 1618  715 (S29W30 L=133)   306     60        0 (1583/ 692-1641/ 746)  11.3 W27-W32 S30-S27
     14  312  647 (S02E57 L= 46)  1285    197        4 ( 261/ 608- 371/ 685) -24.0 E63-E51 S04-N01
     15  510  355 (S23E56 L= 46)   440     71        1 ( 479/ 331- 547/ 383) -59.0 E58-E54 S25-S20
     16  622  331 (S27E49 L= 54)   452     83        0 ( 579/ 313- 665/ 353)  -7.1 E52-E44 S28-S25
  

    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