Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1597 1578 (N27W51 L=314)  1909    331        7 (1472/1468-1711/1662) -51.3 W43-W57 N18-N33
      2 1145 1536 (N30W15 L=277)   894    280        0 (1039/1432-1218/1618) -61.9 W09-W17 N21-N36
      3 1616 1479 (N20W48 L=310)   328     37        0 (1590/1461-1646/1496) -26.6 W45-W50 N19-N21
      4  731 1440 (N28E14 L=248)   651    143        0 ( 680/1396- 804/1484) -37.9 E18-E09 N24-N31
      5 1830  907 (S16W57 L=320)   463     73        0 (1810/ 875-1862/ 932)  31.0 W54-W60 S18-S14
      6 1719  622 (S33W50 L=312)  1366    249        2 (1628/ 586-1819/ 669)   0.0 W40-W63 S35-S30
      7 1372  556 (S34W19 L=281)   652    108        0 (1344/ 506-1393/ 623) -83.7 W16-W20 S37-S28
      8 1060  438 (S38E06 L=257)   337     86        0 (1026/ 410-1085/ 482)  26.6 E08-E02 S39-S34
      9 1287  457 (S40W12 L=275)   437     88        1 (1242/ 437-1335/ 476)   0.0 W08-W15 S41-S38
     10  800  422 (S35E25 L=238)   370     75        0 ( 771/ 395- 835/ 451) -53.1 E26-E23 S37-S33
     11  870  231 (S52E30 L=232)   329     44        0 ( 838/ 217- 899/ 244)   0.0 E33-E26 S53-S50
     12  744  162 (S57E55 L=207)   579     35        0 ( 732/ 135- 756/ 187)  15.5 E66-E48 S59-S54
  

    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