Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1596 1334 (N15W44 L=247)   520    102        1 (1551/1310-1632/1363)   0.0 W39-W46 N12-N16
      2 1138 1325 (N23W12 L=215)   503    123        2 (1073/1301-1203/1347) -29.7 W08-W16 N20-N25
      3  861 1278 (N24E06 L=197)   280     53        3 ( 831/1258- 892/1299)  14.0 E08-E04 N22-N25
      4  951 1248 (N21E01 L=202)   341    103        2 ( 903/1207- 998/1288) -51.3 E03-W01 N17-N23
      5  812 1212 (N21E10 L=193)   333    116        1 ( 795/1161- 837/1273) -90.0 E11-E09 N17-N24
      6 1692 1239 (N07W49 L=252)   524     98        0 (1652/1204-1722/1269) -51.3 W45-W51 N04-N09
      7 1588 1084 (N00W37 L=240)   487     96        3 (1549/1055-1628/1114) -31.0 W34-W39 S01-N02
      8  915  983 (N06E07 L=196)   336     88        2 ( 889/ 931- 933/1017) -90.0 E08-E05 N02-N07
      9 1409  826 (S12W21 L=224)   495    161        0 (1393/ 750-1428/ 910)  78.7 W19-W22 S16-S06
     10  490  612 (S12E42 L=161)   322     66        2 ( 467/ 590- 514/ 647)  53.1 E44-E40 S13-S09
  

    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