Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  536 1750 (N52E44 L=259)   460     57        0 ( 516/1721- 556/1778) -78.7 E46-E41 N49-N54
      2 1729 1531 (N26W60 L=  4)  1117    156        1 (1665/1476-1786/1567) -24.4 W53-W66 N22-N28
      3 1900  943 (S11W66 L=  9)   859    150        4 (1852/ 898-1941/1017)  32.5 W60-W73 S13-S06
      4 1144  823 (S15W06 L=309)   500    214        1 (1073/ 742-1204/ 912) -57.5 W02-W09 S20-S08
      5 1601  840 (S16W37 L=340)   495    156        0 (1546/ 778-1652/ 892)  40.6 W33-W40 S19-S13
      6  577  622 (S23E33 L=270)   289     74        0 ( 552/ 586- 602/ 653)  38.7 E36-E30 S25-S21
      7  233  628 (S20E67 L=236)   448     38        1 ( 212/ 604- 259/ 650) -21.8 E69-E62 S21-S18
      8 1701  610 (S31W50 L=353)   422     86        0 (1662/ 574-1743/ 649)  33.7 W46-W54 S33-S28
      9  950  546 (S31E09 L=294)   255     50        0 ( 927/ 526- 975/ 565) -33.7 E10-E06 S32-S29
  

    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