Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1645 1654 (N27W70 L=200)   567     68        1 (1616/1634-1680/1670) -38.7 W66-W74 N24-N28
      2  864 1564 (N43W02 L=133)   230     40        2 ( 837/1540- 889/1583) -14.0 W00-W04 N41-N44
      3 1051 1538 (N37W15 L=145)   277     69        1 (1021/1513-1089/1562)   0.0 W12-W18 N35-N38
      4  817 1480 (N39E03 L=127)   588    170        0 ( 798/1419- 858/1558)  78.7 E05-W00 N34-N44
      5 1518 1372 (N15W40 L=170)   306     86        1 (1494/1331-1537/1405)  78.7 W38-W41 N12-N17
      6  827 1210 (N22E08 L=122)   987    334        0 ( 761/1085- 931/1338) -74.5 E11-E04 N12-N30
      7  149  901 (N13E70 L= 60)  1061    114        2 ( 118/ 856- 176/ 931) -90.0 E76-E65 N10-N15
      8 1050  897 (S01E01 L=129)   211     40        2 (1018/ 877-1076/ 917)   0.0 E02-W00 S01-N00
      9  921  839 (S01E10 L=120)   239     86        1 ( 873/ 820- 961/ 860)   0.0 E12-E07 S02-N00
     10 1636  740 (S24W34 L=164)  1387    283        2 (1554/ 616-1700/ 848)  42.7 W26-W39 S30-S18
     11  495  711 (S02E40 L= 91)   328     87        0 ( 484/ 668- 508/ 758)  59.0 E41-E37 S04-N01
     12 1162  663 (S17W01 L=131)   559     74        3 (1127/ 626-1199/ 692)  26.6 E01-W03 S19-S14
     13  841  562 (S17E21 L=109)   239    123        0 ( 824/ 513- 869/ 609)  74.1 E22-E18 S20-S13
     14  906  519 (S21E19 L=112)   376     70        0 ( 883/ 487- 924/ 557)  63.4 E20-E16 S22-S18
     15 1045  507 (S24E10 L=120)   289    111        0 (1008/ 468-1070/ 553)  45.0 E13-E07 S26-S21
     16  716  514 (S17E31 L=100)   254     83        1 ( 680/ 492- 761/ 534) -36.9 E32-E28 S19-S15
     17  943  311 (S36E25 L=106)  5051    817       13 ( 669/ 167-1230/ 507)  13.0 E65-W01 S44-S28
     18  384  428 (S18E64 L= 67)   632     83        0 ( 344/ 406- 419/ 453)  12.1 E71-E57 S19-S16
     19 1581  359 (S54W32 L=163)   704     64        1 (1556/ 340-1613/ 387) -32.5 W27-W38 S56-S49
  

    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