Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1449 1601 (N28W42 L=272)   469     97        1 (1404/1581-1503/1620) -26.6 W38-W46 N26-N30
      2  968 1596 (N37W07 L=238)   279     24        0 ( 944/1574- 989/1613) -33.7 W05-W08 N35-N38
      3 1243 1490 (N25W22 L=253)   212     47        0 (1208/1477-1269/1501) -36.9 W20-W24 N23-N26
      4  444 1308 (N27E35 L=196)  1475    315        1 ( 340/1223- 594/1385)  15.9 E43-E22 N23-N29
      5  220 1276 (N29E58 L=173)  1518    153        1 ( 147/1241- 295/1306) -14.6 E71-E48 N25-N31
      6 1542  829 (S19W29 L=260)  1690    579        2 (1382/ 650-1720/1010)  29.7 W16-W45 S28-S11
      7  278  817 (N00E54 L=177)   571    118        1 ( 259/ 749- 300/ 873)  50.2 E57-E50 S03-N03
      8  350  609 (S13E54 L=176)   962    141        2 ( 318/ 535- 378/ 669)  71.6 E57-E51 S17-S08
  

    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