Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1049 1797 (N50W34 L=145)   303     44        1 (1030/1765-1068/1822)  26.6 W29-W37 N48-N52
      2 1393 1771 (N36W56 L=167)   547     78        1 (1359/1745-1433/1792) -15.9 W51-W60 N34-N37
      3 1065 1598 (N36W21 L=132)   276     46        0 (1047/1572-1091/1620)  11.3 W18-W23 N34-N37
      4 1425 1396 (N15W34 L=145)   228     37        0 (1401/1380-1453/1414) -33.7 W32-W35 N13-N15
      5 1028 1371 (N24W10 L=121)   293     56        1 ( 997/1351-1064/1388) -26.6 W07-W11 N22-N25
      6 1833 1332 (S01W64 L=175)  1187    107        2 (1780/1296-1884/1358) -24.4 W59-W72 S03-N02
      7 1274 1145 (N05W17 L=128)   184     91        0 (1233/1119-1312/1171)   0.0 W14-W19 N04-N05
      8  489 1124 (N23E30 L= 81)  1567    279        2 ( 349/1092- 625/1164) -24.2 E41-E21 N19-N28
      9 1459  851 (S16W21 L=132)   819    358        0 (1348/ 735-1555/ 975)  25.2 W11-W28 S19-S10
     10 1486  701 (S26W19 L=130)  1168    551        0 (1383/ 497-1608/ 944)  41.2 W08-W32 S36-S14
     11  768  895 (N04E17 L= 94)   495    170        0 ( 699/ 870- 821/ 941) -15.9 E20-E13 N02-N06
     12  144  856 (N14E70 L= 42)  1483    124        4 (  89/ 816- 202/ 895) -11.3 E81-E61 N11-N17
     13  326  711 (N03E53 L= 58)  3121    521        9 ( 226/ 524- 397/ 886)  29.4 E77-E43 S06-N12
     14 1124  871 (S06W02 L=113)   234    111        0 (1075/ 852-1178/ 884) -45.0 W00-W04 S08-S04
     15 1732  727 (S32W40 L=151)  1929    354        2 (1655/ 615-1857/ 884)  27.6 W32-W54 S38-S26
     16  611  629 (S07E34 L= 77)   473    104        1 ( 577/ 586- 658/ 664)  14.0 E38-E30 S09-S05
     17 1320  357 (S44E01 L=111)   358     64        1 (1298/ 330-1340/ 398) -90.0 E02-W00 S47-S40
     18  947  214 (S43E36 L= 76)  1266    287        2 ( 831/ 169-1063/ 293)   6.3 E48-E21 S46-S40
  

    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