Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  922 1579 (N25E22 L=269)  2102    513        3 ( 821/1338- 984/1728) -47.6 E34-E13 N10-N34
      2  534 1549 (N14E46 L=245)   826    117        4 ( 499/1497- 568/1601) -74.1 E48-E43 N10-N17
      3 1777 1204 (N27W47 L=338)  3644    625        2 (1662/1073-1913/1371) -11.0 W34-W66 N19-N34
      4  307 1040 (S20E48 L=244)   534    129        1 ( 286/ 995- 347/1096) -71.6 E49-E43 S22-S16
      5 1450  972 (N04W26 L=317)   952    366        1 (1281/ 916-1516/1090)  52.4 W16-W28 S02-N11
      6 1918  953 (N17W69 L=  0)  1350    136        4 (1889/ 892-1950/1015) -19.7 W62-W78 N14-N19
      7 1869  933 (N14W61 L=353)   511     58        2 (1850/ 899-1886/ 964)  90.0 W59-W63 N12-N16
      8 1828  727 (N02W64 L=356)  1661    168        1 (1775/ 686-1901/ 790)   0.0 W56-W80 N00-N03
      9  973  664 (S26W07 L=298)   634    254        1 ( 875/ 631-1097/ 713)  45.0 W01-W12 S30-S19
     10  388  563 (S49E39 L=252)  2279    287        3 ( 284/ 495- 508/ 640)  -6.3 E53-E24 S53-S45
     11 1587  541 (S15W51 L=342)  1290    333        1 (1445/ 503-1715/ 577)  30.1 W42-W61 S21-S10
     12 1333  205 (S43W63 L=354)  2522    176        2 (1261/ 167-1409/ 274)  -2.0 W49-W77 S45-S40
  

    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