Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1481 1706 (N30W58 L= 83)  1335    211        6 (1394/1663-1551/1746) -56.3 W52-W64 N25-N34
      2 1295 1653 (N33W39 L= 64)   391     54        1 (1267/1632-1323/1678)  76.0 W36-W41 N30-N34
      3  618 1607 (N52E13 L= 12)   349     39        1 ( 594/1580- 635/1632) -78.7 E15-E11 N49-N54
      4  442 1549 (N53E33 L=352)  1839    269        4 ( 329/1481- 543/1607)   6.3 E48-E20 N50-N55
      5  944 1350 (N25W04 L= 29)   203     59        0 ( 925/1317- 961/1381)  76.0 W02-W05 N23-N27
      6  308 1262 (N36E46 L=339)   346     54        1 ( 290/1236- 325/1293)  45.0 E48-E44 N34-N38
      7  957 1165 (N14E00 L= 25)   453    120        1 ( 913/1130-1012/1209) -63.4 E02-W02 N11-N17
      8 1790  926 (S22W48 L= 73)   582    141        0 (1726/ 903-1834/ 965) -10.3 W42-W53 S24-S19
      9  914  638 (S14E16 L=  8)   749    254        1 ( 850/ 540- 984/ 733)  29.7 E23-E09 S18-S10
     10 1175  614 (S22E02 L= 23)   773    234        3 (1113/ 531-1219/ 714)  38.7 E07-W03 S26-S17
     11 1403  513 (S36W10 L= 35)   761    225        2 (1351/ 431-1449/ 617)  45.0 W04-W15 S41-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