Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  712 1433 (N08E29 L=284)  1147    293        1 ( 676/1293- 748/1569) -60.6 E34-E24 N00-N16
      2 1455 1453 (N31W14 L=327)   269     47        0 (1425/1438-1484/1466)  36.9 W11-W15 N29-N32
      3 1684 1339 (N32W33 L=347)   341     57        0 (1660/1312-1708/1367) -14.0 W31-W35 N30-N33
      4  900 1297 (N05E14 L=300)   199     79        0 ( 865/1279- 944/1312)  26.6 E15-E11 N03-N06
      5  406 1149 (S15E41 L=272)  1802    443        2 ( 238/1100- 550/1176)  13.5 E55-E30 S19-S11
      6 1875  746 (N06W70 L= 23)   804    146        0 (1846/ 730-1909/ 763)   9.5 W64-W76 N04-N07
      7 1522  685 (S10W40 L=353)   423    122        0 (1499/ 620-1547/ 741) -66.8 W38-W41 S13-S06
      8 1225  500 (S29W30 L=343)   686    100        0 (1171/ 464-1267/ 545)  81.9 W27-W31 S32-S25
      9 1687  503 (S13W62 L= 16)   555    109        0 (1659/ 477-1725/ 539)  78.7 W58-W64 S15-S10
     10 1090  463 (S36W23 L=337)   257     36        0 (1072/ 442-1108/ 488) -56.3 W21-W25 S37-S34
     11 1523  347 (S27W59 L= 12)  2120    395        2 (1439/ 293-1598/ 394)   2.9 W49-W69 S29-S24
  

    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