Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1596 1437 (N34W39 L= 33)  1520    285        4 (1540/1313-1686/1533) -71.6 W34-W45 N26-N41
      2 1099 1397 (N25W00 L=355)   541    147        1 (1070/1329-1127/1464)  90.0 E01-W01 N20-N29
      3  829 1269 (N13E15 L=339)   265     70        0 ( 793/1245- 864/1299) -21.8 E17-E12 N12-N15
      4  930 1228 (N12E08 L=346)   275     89        1 ( 880/1206- 965/1251)  31.0 E11-E06 N10-N13
      5  180  862 (S20E66 L=289)  2843    281        5 ( 154/ 736- 210/ 987)  86.2 E72-E60 S27-S12
      6 1678  881 (S00W46 L= 40)   329     94        0 (1654/ 840-1696/ 925)  68.2 W44-W47 S03-N02
      7 1708  715 (S10W53 L= 47)   815    194        1 (1635/ 661-1762/ 783)  45.0 W47-W56 S14-S05
      8  663  681 (S25E20 L=334)   263     53        0 ( 647/ 642- 680/ 708)  78.7 E21-E18 S27-S22
      9 1316  518 (S28W28 L= 22)   686    117        1 (1257/ 486-1364/ 554)  50.2 W24-W30 S30-S24
     10 1472  496 (S27W41 L= 35)   892    143        0 (1416/ 455-1524/ 545)  54.5 W36-W43 S30-S23
     11  742  363 (S48E14 L=340)  2428    409        4 ( 695/ 205- 793/ 512)  71.0 E25-E09 S67-S35
     12 1617  453 (S28W57 L= 52)  2108    299        5 (1482/ 402-1747/ 483)  10.6 W46-W79 S33-S25
  

    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