Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1271 1469 (N27W27 L=184)   456    110        1 (1230/1432-1317/1502)  14.0 W22-W30 N26-N28
      2  475 1264 (N32E33 L=124)  1570    440        2 ( 368/1117- 618/1404) -47.7 E43-E23 N20-N42
      3  199 1220 (N32E59 L= 98)   568    112        0 ( 176/1202- 231/1247) -18.4 E61-E55 N31-N34
      4  216 1117 (N26E57 L=100)   855    118        1 ( 196/1073- 253/1165) -80.5 E59-E53 N22-N28
      5 1675 1108 (S03W43 L=200)   857    139        2 (1619/1079-1732/1143) -18.4 W38-W47 S05-N00
      6 1532 1054 (S02W31 L=188)   552    147        3 (1506/ 977-1571/1115)  45.0 W28-W35 S06-N00
      7  671 1071 (N16E21 L=136)   300    128        0 ( 631/1051- 737/1088)   0.0 E23-E16 N15-N17
      8  281  995 (N18E52 L=105)   319     83        0 ( 267/ 952- 292/1039)  59.0 E53-E50 N15-N20
      9  853  902 (N03E12 L=144)  1067    176        2 ( 788/ 810- 927/ 985) -76.0 E14-E08 S03-N09
     10 1545  833 (S16W28 L=185)   470     73        0 (1529/ 788-1561/ 879)  78.7 W26-W29 S18-S13
     11  839  579 (S16E21 L=136)   775    175        2 ( 766/ 506- 888/ 614) -66.0 E24-E18 S21-S12
     12  687  277 (S35E45 L=112)  1030    203        1 ( 645/ 237- 721/ 321)  28.6 E51-E40 S38-S31
  

    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