Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1496 1546 (N33W43 L=355)   983    221        2 (1396/1504-1591/1589) -29.7 W36-W50 N29-N37
      2 1606 1372 (N21W45 L=357)   815    180        0 (1570/1278-1654/1428)  36.0 W39-W51 N15-N23
      3 1800 1265 (N11W61 L= 13)   561     82        1 (1767/1235-1832/1288) -39.8 W57-W64 N08-N13
      4 1588 1219 (N12W40 L=352)   345     76        0 (1562/1191-1615/1247) -53.1 W37-W41 N09-N13
      5 1770  838 (S14W52 L=  5)   383    160        0 (1756/ 773-1791/ 918)  71.6 W51-W54 S18-S09
      6 1681  656 (S25W46 L=358)   368     71        0 (1643/ 633-1711/ 681)  18.4 W42-W48 S26-S23
      7  807  582 (S21E18 L=294)   228     88        0 ( 764/ 559- 847/ 602)  18.4 E21-E15 S22-S19
      8  889  448 (S31E15 L=297)   691    259        0 ( 848/ 342- 950/ 550)  51.3 E21-E09 S39-S24
      9 1095  504 (S29W00 L=312)   270     96        0 (1077/ 459-1118/ 549)  66.8 E01-W02 S32-S25
     10  786  292 (S44E29 L=283)   609     61        1 ( 755/ 261- 820/ 325) -71.6 E31-E26 S47-S40
  

    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