Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1135 1856 (N52E29 L= 92)   925     81        2 (1107/1819-1168/1895) -30.3 E36-E24 N47-N55
      2 1200 1788 (N48E18 L=103)   472     61        2 (1185/1752-1224/1817) -50.2 E20-E15 N44-N51
      3  803 1741 (N30E40 L= 81)   869    251        2 ( 726/1642- 891/1809) -14.6 E51-E28 N26-N34
      4  332 1378 (S03E56 L= 65)   431     70        0 ( 308/1352- 351/1416)  90.0 E57-E53 S05-N00
      5  122 1116 (S22E75 L= 45)   802     74        1 ( 109/1070- 139/1147)  51.3 E79-E72 S24-S19
      6 1607  791 (S01W41 L=162)  1749    385        2 (1527/ 644-1653/ 964) -32.3 W30-W49 S07-N05
      7  206  864 (S36E60 L= 61)  1400    159        2 ( 167/ 798- 255/ 930)  33.7 E68-E52 S41-S31
      8  742  790 (S27E10 L=111)   326    119        0 ( 688/ 764- 791/ 827)  50.2 E12-E07 S29-S23
      9  549  645 (S41E21 L=100)   305     60        1 ( 517/ 624- 578/ 665)   0.0 E23-E18 S42-S39
     10 1180  517 (S29W27 L=147)   267     83        0 (1144/ 501-1230/ 532)  59.0 W24-W28 S31-S26
     11 1353  275 (S37W57 L=178)  1399    145        3 (1280/ 233-1419/ 313)  77.5 W51-W61 S41-S32
     12  965  157 (S62W51 L=171)   475     45        0 ( 939/ 137- 996/ 180) -45.0 W45-W56 S65-S58
  

    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