Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  735 1865 (N38E57 L=136)   661     51        1 ( 706/1842- 772/1883)   0.0 E63-E52 N35-N39
      2  574 1731 (N24E55 L=138)   386     40        2 ( 555/1708- 594/1756) -71.6 E57-E53 N22-N25
      3  821 1616 (N22E30 L=163)   793    241        3 ( 774/1502- 878/1711) -36.3 E37-E22 N16-N28
      4  407 1519 (N07E55 L=138)  1862    369        4 ( 338/1383- 504/1613)  90.0 E60-E50 S02-N14
      5 1116 1534 (N26E10 L=183)   339    124        0 (1098/1474-1136/1586) -60.3 E11-E06 N22-N29
      6 1210 1306 (N15W03 L=195)   335    115        0 (1176/1250-1225/1349) -51.3 W00-W04 N11-N16
      7  416 1066 (S19E38 L=154)  1687    311        3 ( 291/ 947- 512/1171)  52.6 E48-E33 S27-S10
      8  149 1099 (S22E68 L=125)   864    107        4 ( 109/1060- 192/1142)  -7.6 E75-E60 S23-S19
      9  515  836 (S30E27 L=166)  2726    659        5 ( 231/ 725- 822/ 922)  27.5 E56-E08 S43-S17
     10 1661  837 (N03W43 L=236)   729    107        3 (1622/ 806-1709/ 879)  15.9 W39-W46 N00-N04
     11 1423  792 (S07W29 L=222)   217     93        0 (1396/ 756-1464/ 822)  68.2 W28-W30 S09-S04
     12 1163  617 (S24W20 L=213)   344    104        0 (1109/ 598-1212/ 645)  45.0 W17-W22 S26-S21
     13 1649  551 (S12W56 L=249)   914    195        3 (1596/ 487-1704/ 606)  83.7 W53-W58 S16-S07
     14  568  471 (S53E16 L=176)   409    108        0 ( 535/ 433- 614/ 516) -16.7 E20-E10 S54-S50
     15 1643  438 (S18W64 L=257)   541     57        2 (1622/ 405-1666/ 465) -53.1 W60-W66 S19-S15
     16 1022  293 (S49W31 L=224)  1883    345        5 ( 863/ 223-1143/ 386)   0.0 W12-W46 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