Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1233 1449 (N30W10 L=287)   219     50        0 (1218/1417-1247/1478) -76.0 W08-W11 N28-N32
      2 1443 1245 (N19W25 L=303)   568    185        0 (1371/1180-1517/1277) -21.8 W20-W30 N15-N20
      3 1800 1117 (N14W56 L=333)   715    176        0 (1763/1045-1830/1197) -60.9 W52-W59 N09-N18
      4  191 1005 (S09E63 L=214)   490     75        0 ( 155/ 986- 231/1022)  -5.2 E68-E57 S09-S07
      5  165  909 (S15E67 L=210)  1017    129        0 ( 117/ 859- 215/ 942) -14.9 E77-E60 S17-S13
      6  464  860 (S14E36 L=242)   323     65        1 ( 431/ 839- 491/ 890)  21.8 E38-E33 S15-S12
      7 1037  850 (S09W02 L=280)   289     69        1 (1009/ 815-1065/ 882)  63.4 W00-W04 S11-S07
      8  288  833 (S18E52 L=225)   369     75        0 ( 271/ 793- 305/ 863) -90.0 E53-E49 S20-S16
      9  329  754 (S23E49 L=229)   800    116        2 ( 288/ 713- 386/ 792)  -5.7 E53-E43 S24-S20
     10  290  708 (S27E54 L=223)   393     46        1 ( 267/ 686- 316/ 727) -33.7 E56-E50 S27-S25
     11 1705  570 (S21W59 L=336)   412     66        0 (1680/ 540-1725/ 599)  63.4 W57-W60 S22-S18
  

    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