Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  473 1633 (N15E55 L=278)  1014    188        1 ( 421/1557- 521/1713) -85.2 E56-E52 N09-N21
      2 1154 1253 (N09W01 L=333)   470    143        1 (1104/1182-1197/1303) -35.5 E02-W05 N05-N10
      3  498 1208 (S09E36 L=296)   322     98        0 ( 461/1170- 525/1235) -36.9 E38-E33 S11-S08
      4  420 1194 (S12E41 L=291)   320     66        1 ( 393/1173- 452/1215)  53.1 E43-E39 S13-S09
      5 1547 1092 (N11W28 L=  0)   433    115        0 (1531/1030-1561/1154) -56.3 W26-W30 N07-N14
      6 1396  993 (N01W21 L=354)   706    253        0 (1316/ 877-1434/1091) -30.3 W14-W26 S04-N04
      7  695  811 (S27E14 L=318)   697    289        2 ( 566/ 782- 840/ 837)  14.7 E23-E04 S30-S24
      8  412  758 (S38E35 L=298)   280     82        1 ( 376/ 746- 451/ 774)  18.4 E37-E31 S38-S36
      9 1624  732 (S05W45 L= 17)   346     82        0 (1592/ 696-1647/ 760)  -9.5 W41-W47 S06-S04
     10  889  551 (S38W06 L=338)   256     74        0 ( 853/ 527- 922/ 582)  59.0 W04-W08 S39-S34
     11 1548  384 (S25W57 L= 29)   462    102        0 (1519/ 369-1585/ 407)   8.1 W53-W60 S26-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