Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  881 1764 (N60W08 L= 14)   501     94        0 ( 851/1733- 915/1786)   0.0 W04-W13 N56-N62
      2  801 1732 (N59E02 L=  3)   739     99        1 ( 753/1698- 833/1756) -36.9 E08-W01 N55-N61
      3 1234 1717 (N46W37 L= 43)  1186    178        1 (1147/1689-1320/1747)  -6.3 W28-W46 N43-N48
      4  726 1696 (N58E12 L=354)   817     87        3 ( 694/1648- 752/1732) -82.9 E15-E08 N53-N61
      5 1463 1688 (N38W55 L= 61)   786    107        1 (1421/1663-1519/1704) -10.3 W49-W60 N35-N39
      6 1617 1446 (N19W51 L= 57)   306     99        0 (1587/1402-1642/1476) -78.7 W49-W52 N16-N21
      7 1367 1261 (N14W26 L= 32)  1353    267        3 (1307/1153-1406/1366)  81.3 W21-W28 N08-N21
      8  639 1339 (N32E21 L=345)   283     78        0 ( 609/1309- 681/1359)  18.4 E23-E17 N30-N33
      9  624 1089 (N17E24 L=342)   441    142        0 ( 561/1062- 688/1117)   6.3 E28-E19 N15-N17
     10 1818  703 (S31W57 L= 63)  2065    166        3 (1782/ 643-1853/ 782) -83.7 W52-W63 S34-S25
     11  502  646 (S08E41 L=325)   275     77        0 ( 463/ 630- 542/ 662) -21.8 E43-E38 S09-S07
     12  992  239 (S47E23 L=343)   918    156        0 ( 926/ 203-1063/ 284)  11.3 E29-E14 S49-S43
  

    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