Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1274 1681 (N36W35 L=177)  1332    269        1 (1179/1616-1384/1715) -14.0 W26-W43 N32-N37
      2  568 1544 (N45E21 L=121)   337    112        1 ( 515/1522- 624/1562)  -5.7 E25-E15 N43-N45
      3 1479 1043 (S06W27 L=169)  1081    355        2 (1403/ 872-1550/1211)  43.2 W19-W35 S14-N01
      4 1405  734 (S22W17 L=159)   740    276        1 (1275/ 696-1518/ 784) -17.4 W08-W24 S24-S17
      5 1066  471 (S31E11 L=131)   969    233        2 ( 991/ 394-1098/ 567)  54.0 E17-E07 S35-S24
      6  651  426 (S25E39 L=103)   857    325        1 ( 628/ 316- 667/ 558)  53.1 E45-E32 S32-S16
      7 1495  526 (S39W21 L=163)   384     57        2 (1462/ 500-1524/ 550)  21.8 W17-W23 S40-S36
  

    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