Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  583 1641 (N38E44 L=328)   367     54        2 ( 559/1612- 604/1665) -23.2 E48-E41 N36-N40
      2  739 1621 (N39E29 L=343)  1048    155        1 ( 672/1584- 815/1657)  -8.1 E35-E21 N36-N40
      3  315 1521 (N27E67 L=305)   515     42        1 ( 291/1503- 345/1534)   0.0 E72-E62 N25-N27
      4 1316 1392 (N27W17 L= 30)  1312    346        1 (1242/1294-1371/1499) -58.4 W11-W21 N20-N34
      5  940 1463 (N29E10 L=  3)   213     41        1 ( 919/1448- 964/1485) -33.7 E10-E07 N28-N30
      6  471 1366 (N19E42 L=330)   582    181        0 ( 405/1298- 512/1400) -22.6 E48-E36 N14-N20
      7 1754 1267 (N22W54 L= 66)   432     99        0 (1731/1214-1771/1311)  80.5 W51-W55 N18-N24
      8 1573 1194 (N16W36 L= 48)   241     63        0 (1555/1160-1592/1215)  63.4 W34-W37 N13-N17
      9  937 1101 (N06E06 L=  6)   579    152        1 ( 883/1053- 979/1145)  50.2 E08-E03 N02-N08
     10 1317  562 (S25W24 L= 36)   440    109        0 (1260/ 542-1368/ 584)   0.0 W19-W27 S26-S23
  

    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