Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  405 1650 (N24E68 L=230)  1355    108        3 ( 377/1597- 438/1694) -25.0 E81-E61 N19-N26
      2  617 1643 (N27E45 L=253)   382    102        1 ( 574/1606- 654/1664)  -6.3 E49-E40 N25-N28
      3 1391 1434 (N31W15 L=313)   406    155        0 (1340/1372-1444/1494) -29.1 W10-W19 N28-N34
      4 1225 1437 (N27W04 L=302)   870    328        1 (1099/1336-1342/1492) -12.5 E05-W13 N23-N29
      5 1888 1183 (N27W67 L=  5)  2099    243        0 (1850/1078-1908/1285) -81.3 W62-W71 N20-N33
      6  471 1059 (S11E36 L=263)   419     79        1 ( 427/1039- 506/1080)  31.0 E38-E33 S12-S09
      7 1785  986 (N12W53 L=351)   736    124        0 (1759/ 930-1807/1041)  76.0 W51-W55 N08-N16
      8 1467  537 (S21W42 L=340)  2199    405        2 (1388/ 421-1574/ 637)  82.9 W37-W51 S28-S12
      9  487  544 (S44E31 L=267)   362     37        1 ( 473/ 514- 507/ 575)  68.2 E32-E28 S46-S41
     10  930  519 (S35W06 L=304)   593    175        0 ( 848/ 491-1014/ 549)  26.6 W00-W11 S37-S30
  

    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