Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1414 1824 (N48W71 L=327)   794     72        0 (1386/1797-1441/1851)   0.0 W62-W88 N47-N49
      2 1533 1682 (N36W61 L=317)  1207    236        0 (1491/1613-1578/1732)  56.3 W56-W69 N30-N39
      3  717 1517 (N43E14 L=242)   508    124        2 ( 668/1493- 793/1543) -20.0 E18-E07 N40-N45
      4 1564 1446 (N21W46 L=302)   321     83        0 (1540/1413-1594/1468)   9.5 W43-W49 N19-N22
      5 1752 1380 (N13W61 L=317)   574    121        2 (1721/1324-1784/1434) -90.0 W59-W63 N08-N16
      6  954 1370 (N29W02 L=257)   202     85        0 ( 913/1354- 992/1387) -36.9 W00-W04 N27-N30
      7 1270 1288 (N18W20 L=276)  1133    397        0 (1089/1230-1429/1333) -31.0 W09-W29 N11-N23
      8  476 1307 (N32E35 L=221)   316    123        1 ( 423/1281- 528/1326)  -6.3 E39-E30 N30-N33
      9  217 1140 (N24E61 L=195)  1851    234        4 ( 186/1061- 259/1239)  73.3 E65-E55 N18-N29
     10  365 1130 (N22E45 L=211)   285     89        0 ( 334/1117- 419/1151) -20.6 E47-E39 N20-N23
     11 1727 1008 (S08W47 L=303)   511     83        1 (1705/ 963-1745/1045)  78.7 W44-W48 S10-S05
     12  167  918 (N10E69 L=187)   620     80        2 ( 140/ 877- 183/ 943) -53.1 E73-E66 N07-N11
     13  139  875 (N07E75 L=181)  1029     67        1 ( 116/ 846- 158/ 907) -26.6 E80-E71 N05-N08
     14 1595  415 (S47W38 L=294)   653     82        1 (1557/ 390-1627/ 446)  -7.1 W32-W42 S49-S44
     15  592  317 (S34E51 L=205)   929     28        1 ( 551/ 271- 626/ 381)  27.9 E61-E44 S37-S28
     16 1299  256 (S54W07 L=262)   484     81        1 (1274/ 233-1328/ 281)  21.8 W03-W09 S56-S51
  

    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