Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  903 1371 (N28E05 L=158)   229     58        1 ( 874/1349- 935/1386)   0.0 E07-E03 N26-N28
      2 1537 1283 (N16W37 L=200)   519    156        1 (1465/1263-1608/1305) -21.8 W31-W42 N14-N18
      3  689 1044 (N08E21 L=142)   413    120        0 ( 654/1017- 741/1091) -38.7 E23-E17 N06-N11
      4 1894 1013 (S05W68 L=231)  1313     83        2 (1857/ 988-1943/1040) -13.2 W62-W80 S07-S02
      5 1800  825 (S15W56 L=219)   611    122        1 (1758/ 784-1833/ 850) -16.7 W51-W61 S18-S13
      6 1429  702 (S19W24 L=187)   753    127        2 (1372/ 666-1493/ 739)   7.1 W20-W28 S20-S16
      7  678  697 (S13E25 L=138)   405    157        0 ( 618/ 677- 728/ 725) -15.9 E28-E21 S14-S11
      8 1742  638 (S27W53 L=216)   601    116        0 (1706/ 595-1782/ 686)  45.0 W50-W57 S29-S24
      9 1256  527 (S29W12 L=175)  1036    304        2 (1172/ 435-1348/ 645)  47.3 W05-W18 S35-S21
     10  584  327 (S40E47 L=116)   473     52        2 ( 567/ 296- 600/ 365)  36.9 E51-E42 S42-S36
     11 1054  266 (S49E06 L=157)   527     68        0 (1025/ 241-1087/ 295) -31.0 E08-E02 S51-S45
  

    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