Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  748 1864 (N38E56 L=136)   771     53        1 ( 719/1837- 786/1884)  -6.3 E62-E51 N35-N39
      2  685 1814 (N32E54 L=138)   614     52        2 ( 665/1781- 709/1846) -45.0 E58-E50 N29-N35
      3  835 1609 (N22E29 L=163)   775    251        3 ( 790/1494- 896/1712) -38.7 E36-E21 N16-N28
      4  419 1536 (N08E54 L=137)  2956    345        7 ( 345/1404- 530/1623)  75.1 E61-E49 N00-N15
      5 1131 1527 (N26E09 L=183)   379    116        1 (1114/1466-1151/1582) -54.5 E10-E05 N22-N29
      6 1225 1295 (N14W04 L=195)   323    123        0 (1187/1239-1239/1344) -45.0 W00-W05 N11-N16
      7  763 1302 (N02E22 L=170)   200     67        1 ( 743/1270- 782/1336) -76.0 E22-E21 N00-N04
      8  431 1063 (S19E37 L=154)  1645    323        4 ( 302/ 944- 530/1163)  50.5 E46-E31 S27-S10
      9  137 1116 (S21E70 L=122)   655    103        3 ( 111/1078- 171/1156)  29.1 E75-E64 S23-S18
     10  546  836 (S29E24 L=167)  2711    674        3 ( 254/ 724- 843/ 917)  29.9 E53-E06 S43-S16
     11 1636  834 (N02W42 L=233)   270     68        0 (1623/ 803-1654/ 872) -45.0 W39-W43 N00-N03
     12 1515  754 (S06W36 L=228)   295     94        0 (1493/ 706-1535/ 794) -38.7 W33-W38 S08-S04
     13 1185  613 (S24W22 L=214)   374    124        1 (1124/ 595-1248/ 635)  50.2 W18-W25 S26-S20
     14 1693  568 (S10W58 L=250)   415     69        1 (1678/ 524-1712/ 598) -45.0 W55-W61 S12-S08
     15 1633  528 (S14W56 L=248)   706     90        1 (1604/ 482-1664/ 578) -71.6 W53-W58 S16-S10
     16  574  468 (S53E16 L=176)   404    120        0 ( 543/ 428- 620/ 519) -24.0 E19-E10 S54-S49
     17 1649  434 (S18W65 L=257)   527     68        0 (1627/ 398-1672/ 456) -63.4 W61-W67 S20-S16
     18 1029  288 (S49W32 L=223)  1915    357        5 ( 870/ 219-1150/ 382)   0.0 W13-W47 S50-S45
  

    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