Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1757 1360 (N05W59 L=153)   374     81        1 (1746/1323-1772/1399)  45.0 W56-W61 N03-N07
      2  340 1228 (N34E43 L= 50)   518    123        3 ( 318/1185- 373/1263)  53.1 E45-E40 N31-N36
      3  442 1002 (N18E37 L= 56)   543    200        2 ( 375/ 936- 517/1047) -45.0 E41-E32 N13-N22
      4  322  868 (N13E50 L= 44)   324     87        0 ( 287/ 841- 351/ 897) -38.7 E52-E47 N10-N14
      5 1236  726 (S15W05 L= 98)   846    276        3 (1156/ 637-1281/ 792)  53.1 W00-W07 S20-S10
      6 1757  723 (S32W43 L=136)   754    132        3 (1728/ 663-1781/ 787)  60.3 W39-W45 S35-S27
      7  359  617 (S03E55 L= 38)  1157    329        4 ( 303/ 495- 417/ 753)  66.4 E60-E49 S11-N05
      8 1162  707 (S14W00 L= 94)   251     62        2 (1136/ 687-1187/ 731)  18.4 E01-W02 S15-S12
      9 1383  635 (S25W12 L=106)   262    134        0 (1331/ 616-1434/ 679) -45.0 W09-W15 S26-S20
     10  832  254 (S37E38 L= 55)   539     81        3 ( 798/ 235- 868/ 276)   0.0 E40-E34 S38-S35
  

    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