Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  669 1670 (N25E43 L=222)  2010    386        6 ( 543/1563- 807/1741)  64.8 E47-E36 N16-N33
      2  839 1666 (N29E31 L=234)   552    101        2 ( 806/1619- 871/1705) -32.0 E35-E27 N26-N31
      3 1296 1428 (N26W06 L=271)  1979    502        4 (1186/1191-1357/1599) -50.7 E05-W13 N12-N35
      4 1209 1198 (N10W06 L=271)   203     97        1 (1176/1163-1238/1239) -31.0 W03-W08 N08-N11
      5 1616 1048 (N13W35 L=300)   325     99        0 (1602/1001-1627/1108) -68.2 W32-W36 N10-N15
      6  263 1039 (S21E52 L=213)   864    190        2 ( 212/ 958- 329/1086)  48.0 E57-E46 S26-S16
      7  824  987 (S11E10 L=255)   280    182        1 ( 792/ 904- 868/1035)  83.7 E11-E08 S16-S07
      8  641  922 (S19E20 L=245)   382    194        0 ( 612/ 866- 683/ 983) -45.0 E22-E16 S21-S14
      9  310  928 (S26E46 L=219)   670    150        1 ( 236/ 904- 372/ 964)  23.2 E54-E40 S29-S22
     10 1775  837 (N06W54 L=319)  1005    222        2 (1712/ 759-1815/ 937)  79.7 W49-W56 N00-N11
     11  603  459 (S50E15 L=250)  1545    276        3 ( 505/ 390- 716/ 549)  -9.9 E25-E02 S53-S46
  

    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