Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1176 1800 (N51W44 L= 62)   435     73        0 (1144/1772-1214/1820)   5.2 W39-W50 N50-N52
      2  232 1265 (N36E60 L=318)   697     79        1 ( 203/1226- 263/1299) -24.4 E66-E54 N33-N38
      3  772 1085 (N16E14 L=  4)   944    241        3 ( 735/ 998- 825/1199)  85.6 E15-E11 N09-N22
      4  349 1074 (N22E45 L=333)   492    149        0 ( 284/1045- 409/1123) -28.6 E51-E39 N20-N26
      5  894  944 (N05E09 L=  9)   412    117        1 ( 830/ 926- 947/ 972) -39.8 E12-E06 N03-N08
      6  403  897 (N11E43 L=335)   530    129        2 ( 365/ 856- 448/ 935) -38.7 E45-E39 N08-N13
      7 1515  844 (S15W26 L= 44)   235     87        0 (1490/ 797-1533/ 878)  36.9 W23-W27 S17-S12
      8  198  813 (N08E66 L=312)   657     67        0 ( 169/ 791- 233/ 840) -20.6 E69-E61 N06-N10
      9 1087  571 (S21E06 L= 11)   463    163        0 (1049/ 508-1126/ 643)  41.2 E10-E02 S24-S17
     10  574  286 (S31E58 L=320)  1379    184        3 ( 526/ 215- 611/ 372)  51.3 E66-E52 S38-S23
     11  726  329 (S30E40 L=338)   418     77        0 ( 705/ 293- 754/ 352)  33.7 E42-E36 S33-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