Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1673 1606 (N22W65 L=129)  2532    356        5 (1581/1557-1767/1656) -36.3 W58-W75 N16-N27
      2 1466 1383 (N17W37 L=101)   540    194        2 (1413/1340-1494/1459)  29.1 W31-W40 N14-N21
      3  491 1184 (N27E31 L= 33)   361     90        0 ( 464/1150- 525/1217) -51.3 E33-E29 N24-N29
      4 1137  708 (S14E00 L= 64)   683    166        2 (1101/ 627-1176/ 789)  66.0 E03-W03 S18-S09
      5  311  610 (S05E60 L=  4)   786    114        1 ( 288/ 558- 343/ 651)  71.6 E62-E56 S07-S01
      6  461  583 (S08E47 L= 17)  1327    253        3 ( 368/ 542- 551/ 611) -24.8 E55-E40 S12-S05
      7  460  365 (S23E59 L=  5)  2682    420        6 ( 425/ 284- 495/ 455)  51.8 E67-E54 S30-S16
      8  797  333 (S31E34 L= 30)   286    108        0 ( 768/ 305- 818/ 371)  29.7 E37-E30 S32-S28
  

    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