Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  548 1662 (N25E53 L=327)  1001    149        0 ( 477/1626- 602/1687)   9.5 E58-E46 N21-N26
      2  392 1509 (N12E57 L=323)   645    101        1 ( 358/1469- 429/1553)  81.9 E58-E54 N08-N15
      3 1290 1380 (N24W08 L= 28)   304     97        0 (1240/1357-1332/1398)   9.5 W04-W10 N23-N25
      4 1428 1123 (N12W22 L= 42)   204     49        0 (1408/1093-1444/1151) -45.0 W19-W23 N10-N13
      5 1851  960 (N14W60 L= 81)   914    149        2 (1807/ 886-1879/1014) -45.0 W54-W65 N09-N16
      6 1146  440 (S35W25 L= 45)   727    146        1 (1084/ 406-1210/ 479)  -4.8 W19-W31 S35-S32
      7 1560  426 (S24W57 L= 77)   571     63        0 (1522/ 409-1599/ 442)   5.7 W52-W62 S25-S22
      8 1239  251 (S45W47 L= 67)   354     47        0 (1222/ 234-1262/ 276)   0.0 W44-W50 S46-S43
  

    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