Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1400 1816 (N47W50 L=358)   717     85        2 (1363/1785-1425/1842)  12.1 W43-W57 N44-N49
      2 1106 1617 (N35W13 L=322)   311     74        0 (1073/1602-1143/1632) -36.9 W11-W15 N33-N36
      3  448 1438 (N31E36 L=272)   713    150        2 ( 401/1377- 502/1506) -42.3 E42-E31 N26-N36
      4 1676 1367 (N13W48 L=356)   946    156        2 (1619/1315-1741/1427) -36.9 W44-W52 N09-N16
      5  707 1395 (N25E16 L=292)   242     73        0 ( 674/1372- 747/1417)  21.8 E18-E13 N23-N25
      6 1862 1024 (S09W60 L=  8)  2761    295        6 (1779/ 926-1920/1156) -70.3 W51-W67 S14-N00
      7  512 1050 (N06E31 L=277)  1799    396        7 ( 336/ 996- 694/1098)  -4.6 E44-E19 N02-N08
      8 1536  691 (S26W31 L=339)  4766    859        6 (1135/ 582-1871/ 839)   7.7 W02-W62 S30-S20
      9 1171  770 (S17W06 L=315)   187    101        0 (1122/ 742-1218/ 784)   8.1 W02-W09 S18-S16
     10  563  721 (S14E32 L=276)   308    114        1 ( 509/ 702- 621/ 738) -14.0 E36-E28 S14-S12
     11  659  703 (S16E26 L=282)   201     87        0 ( 631/ 660- 676/ 733) -78.7 E27-E25 S18-S13
     12  927  628 (S23E11 L=298)   374    112        1 ( 896/ 595- 956/ 675)  51.3 E12-E08 S25-S20
     13  287  597 (S18E60 L=248)  1298    114        4 ( 249/ 537- 334/ 645)  45.0 E66-E55 S21-S14
     14 1717  500 (S42W55 L=  3)  1501    230        4 (1608/ 455-1819/ 544)   0.0 W41-W73 S43-S39
     15  795  373 (S39E27 L=282)   741    290        0 ( 695/ 323- 909/ 466)  20.0 E37-E15 S41-S33
     16  416  413 (S31E58 L=250)  1666    185        3 ( 354/ 359- 506/ 440)   0.0 E71-E49 S34-S28
  

    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