Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  783 1848 (N37E49 L=343)  1168    207        3 ( 726/1761- 841/1920) -48.8 E58-E43 N29-N45
      2  544 1690 (N20E54 L=339)   548     89        4 ( 503/1669- 583/1719)  26.6 E57-E50 N17-N22
      3  761 1606 (N19E33 L=359)   500    145        2 ( 688/1566- 809/1636)  -5.7 E38-E28 N17-N20
      4  506 1476 (N06E44 L=348)  1218    326        5 ( 358/1431- 621/1566)  40.6 E55-E36 N00-N13
      5 1735 1251 (N27W39 L= 71)   755    142        0 (1686/1202-1783/1284)   0.0 W34-W44 N25-N29
      6  230 1198 (S15E59 L=333)   963    139        2 ( 181/1150- 284/1245)  33.7 E63-E53 S18-S11
      7 1077 1194 (N04E02 L= 30)   333    101        1 (1029/1178-1123/1213)  26.6 E03-W01 N02-N05
      8  201 1029 (S26E59 L=333)   810    194        0 ( 180/ 943- 238/1115) -51.3 E62-E53 S30-S20
      9  389  884 (S30E38 L=354)  1200    217        6 ( 325/ 845- 485/ 939)   3.8 E44-E29 S33-S27
     10 1313  862 (S07W20 L= 52)   203     57        0 (1293/ 833-1327/ 893)  90.0 W19-W21 S09-S05
     11  824  855 (S21E07 L= 25)   244     84        0 ( 785/ 842- 874/ 867)  36.9 E09-E03 S22-S19
     12 1230  714 (S18W21 L= 53)   361     95        1 (1187/ 693-1258/ 757)  68.2 W18-W22 S19-S14
  

    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