Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1175 1733 (N52W22 L= 51)   329     60        0 (1143/1718-1214/1744)   0.0 W17-W25 N51-N53
      2 1563 1562 (N35W50 L= 79)  1714    318        2 (1418/1524-1689/1627) -20.6 W37-W63 N30-N41
      3  546 1425 (N32E34 L=355)   910    218        2 ( 494/1367- 618/1508) -32.7 E41-E27 N28-N37
      4 1356 1303 (N20W24 L= 53)   236    117        0 (1335/1268-1382/1353) -90.0 W22-W26 N17-N22
      5  542 1202 (N17E32 L=358)   386    149        1 ( 504/1142- 596/1244)  45.0 E34-E27 N13-N20
      6  702 1132 (N13E20 L=  9)   591    260        1 ( 615/1032- 774/1217) -50.2 E26-E15 N06-N18
      7 1048 1104 (N10W02 L= 31)   445     98        0 (1013/1063-1078/1147) -59.0 W00-W03 N06-N12
      8 1396  786 (S12W23 L= 52)   431    132        1 (1353/ 750-1458/ 825)  15.9 W19-W26 S14-S10
      9  399  735 (S11E45 L=344)   981    191        2 ( 365/ 653- 464/ 779)  32.5 E50-E39 S16-S08
  

    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